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I.  IN  IBCCOCTION 


GC1SAVG  is  a  simplified  laser  propagation  program.  It  is 
intended  specifically  ior  near  vertical  ground-to-space 
applications  utilizing  a  fixed  earth-based  transmitter 
directed  at  a  sirgle  target  satellite.  The  primary  purpose 
cf  the  program  is  tc  provide  irradiance  and  fluence  cn 
target  alcrg  with  related  propagation  data. 

lie  program  was  written  by  Charles  E.  Hogge  from  the  Air 
force  Reapers  laboratory,  Kirtland  Air  Force  Base, 
Alnuguergue,  New  Mexico.  GU1SAVG  is  one  of  a  family  cf 
ground- tc-space  propagation  pregrams.  Other  versions  include 
GUTSff  (footprint)  ard  GUTSMTF.  GUTSFP  computes  an  engage¬ 
ment  envelope  based  cn  user  supplied  irradiance  threshold 
levels.  She  propagation  calculation  methods  are  identical  to 
GUISA^G.  GO  ISM IF  is  a  full  wave-optics  program  using  fast 
Fourier  transfoms  ir  the  beam  propagation  computa  tiens . 

Ihe  tasic  approach  used  in  GUTSAVG  is  to  utilize  modula¬ 
tion  transfer  functicxs  to  characterize  effects  such  as  team 
guality,  jitter,  turbulence  and  to  apply  these  effects  at 
the  aperture  as  a  sirgle  phase  screen.  Ihe  same  approach  is 
used  in  ESI-IV  [Bef.  1].  The -general  modulation  transfer 
function,  (MTF)  ,  usee  in  this  way  is  described  in  following 
sections  along  with  the  development  of  each  beam  degrading 
mechanism.  A  linearized  model  is  used  for  thermal  blooming 
which  is  also  effectively  applied  at  the  aperture  as  a  phase 
screen.  In  the  case  cf  blooming,  however,  the  phase  variance 
due  tc  blooming  is  computed,  and  the  Strehl  relation  is  used 
to  determine  the  relative  irradiance  reduction.  The  effect 
cf  thermal  blooming  and  the  other  spreading  effects  are 
combined  using  both  IBS  and  multiplicative  methods,  finally, 
the  average  of  these  two  methods  is  used  to  determine  the 


1C 


total  s^stea  irradiance  reducticc  from  tne  dittracticr. 


limited  case. 


II 


PROPAGATION  FEATURES 


A.  TEE  Elf  APPROACH 

G11SA7G  uses  HTFs  to  apply  and  tnen  to  consolidate  tie 
effects  cf  fitter,  beau  quality,  and  turbulence. 

Essentially,  these  effects  are  replaced  by  a  phase  screen  at 
the  aperture  which  multiplies  the  initial  complex  aperture 
distribution  by  a  raiccm  phase  distortion  factor. 

lute nirsxi  and  Ytra  £Ref.  2]  have  developed  an  expres¬ 
sion  fcr  the  average  intensity  at  a  pornt  P  in  the  far-field 
due  tc  perturbations  at  the  aperture.  Ine  intensity  at  P  is 


I  (p)  = 


kY  fH*u  ■z>e*p[-(it)pp]'  dp 

j  U(r+V 


iP)  U*  (r-Jgp)  exp  - 


i  •  M  rp 


2 

dr 


(2.1) 


where 


P  =  |r2  -  r i 


(2.2) 


and 


_  =  r2  +  r ! 
2 


(2.3) 


p"  is  the  vector  from  the  z  axis  of  symmetry  in  the  far-field 
to  the  pcint  P.  U{r)  is  the  aperture  distribution  and  ^("P 
,z)  is  the  MTF  fcr  the  disturbance.  r,  and  r2  are  two  arti- 
trari  pcints  in  the  aperture  where  the  phase  perturbation  is 
ireasuied.  "p  is  the  distance  between  the  points.  z,  in  the 
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case  considered  here,  is  a  constant.  If  the  phase  distur¬ 
bance  is  a  random  Gaussian  variable  with  a  known  correlation 
function,  then  the  Mlf  can  he  expressed  as 

M^(p)  =  <  exp  [i  (r  i  -r  2  )]>  (2.4) 


In  terns  of  the  structure  function. 


M^(p)  =  exp  [(-Jj)D^(r  i-r  2 )] 


Ihe  second  integral  in  equation  2.1  represents  the  urrcrmal- 
ized  aperture  MTI.  Normalizing  this  term  witn  the  power  n 
the  apertuxe,  (P0) ,  tc  cause  the  nlF  to  be  unity  at  the 
origin  results  in 


I  (p) 


fM  (p)Ma  (p)e 


(2.6) 


Kotirc  the  symmetry  of  the  intensity  for  a  given  p  and 
expressing  H  (P)  as  the  combined  effects  of  jitter,  fear 
quality,  and  turbulence,  equation  2.6  can  be  rewritten  as 


I  (p) 


(p)Mt(p)Mb(p)  Ma(p)  J, 


(^f)pdp 


(2.7) 


A  Fcuiier-Eessel  transform  has  teen  used. Note  that  the  tips 
cf  jitter  (hj)  ,  team  qualit  y  (Mj-j  ,  and  turb  ulence  (Mt)  have 
teen  substituted  for 

3he  ESF-IV  manual  £fief.  3]  contains  the  development 
above  in  mere  detail-  MTFs  for  the  specific  effects  are 
described  ir  following  sections. 
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E.  1EEEEAI  BLCOiJING 

lie  thermal  blooming  model  in  GUISAVG  is  based  or.  tie 
following  linearized  density  perturbation  eguation. 


x ' 


2.8 


-co 


X 


Here,  a  { £)  is  the  atmcspher ic  absorption  coefficient  at  a 
distarce  £  along  the  beam,  V 0  is  a  constant  transverse  kind 
velocity  and  p 0  is  the  ambient  pressure.  Ihe  exponential 
term  represents  the  tctal  extinction  due  to  scattering  anc 
absorption,  x'is  a  constant  cf  integration.  It  is  assumed 
that  the  beam  is  propagated  in  the  positive  z  direction  a 
distarce  £  and  that  the  wind  vector  is  in  the  positive  x 
direction.  Ihe  intensity  integral  represents  the  heatinc  cf 
the  atmosphere  as  it  transits  the  beam  £Eef.  4],  Seme 
assumptions  embodied  in  the  above  eguation  are  that  V0<<  c0 
,  the  local  sonic  velocity,  sc  that  the  process  represented 
occurs  at  constant  pressure,  and  that  the  kinetics  cf 
absorption  and  ccnveisicn  to  heat  are  extremely  fast 


£Bef.  5] 


Ey  applying  the  Gladstone-Eale  relation, the  density 
relation  can  be  expressed  as  a  change  in  the  refractive 
index . 


An  =  (n  0-l) 

P  o 


(2.9) 
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Eguaticn  2.6  can  be  rewritten  as 


(2.  10) 


oo 


X  exp 


Ihe  chance  in  the  wavefront  phase  of  a  team  due  to  the 
refractive  index  charge  wner  the  beam  is  propagated  a 
distance  %  •  is  given  iy 


Substituting  equation  2.10  intc  2.11  and  changing  the  limits 
cf  integration  tc  reflect  the  ground- to-space  propagation 

path  results  in  x' 


I (x' , y ) dx ' 


30  km,  ana  hfc  is  the  height  of  the  laser  transmitter.  Also, 
the  wind  term  has  beer  expanded  to  include  the  relative  wind 
velocity  due  to  slewing.  V  0  is  assumed  to  be  parallel  and 
opposite  in  directicr  of  that  cf  the  target  motion.  £  is  the 
angle  cf  incidence  of  the  wind  to  the  beam  so  that  V  0c cs(£) 
represents  the  transverse  wind,  w h  is  the  effective  wind 
generated  hy  slewing.  oj  is  the  angular  slew  rate. 

lie  first  portion  of  eguation  2.12  is  independent  cl 
path  while  the  second  part  is  not,  assuming  I  (x, y)  dees  net 
change  alcnc  the  propagation  path.  This  assumption  .is  valid 
only  fer  a  very  small  amount  of  blooming.  The  approach  used 
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ir.  GUTSAVG  is  to  determine  the  phase  distortion  1  ue  to 
thermal  blooming  by  first  evaluating  the  path  invariant  part 
of  equation  2-  1 2  '.  This  is  accomplished  by  constructing  a 
phase  screen  at  the  aperture  and  then  removing  the  best  fit 
tilt,  focus  curvature,  and  mean  phase.  Zernike  polynomials 
are  used  to  model  these  aberrations.  The  result  is  the  resi¬ 
dual  phase  due  to  thermal  blooming  alone.  The  variance  of 
the  phase  is  then  computed. 

The  path  dependent  term  is  evaluated  within  the  angle 
interval  loop  of  the  program  and  is  applied  to  the  previ¬ 
ously  computed  phase  variance  during  each  path  iteration. 

The  path  iteration  process  is  diagrammed  in  the  engagement 
geometry  section.  Also,  see  Figures  2.1  and  2.2  for  a  flow 
diagram  of  the  general  treatment  of  thermal  blooming  ir  tne 
program. 

Cnee  the  total  phase  variance  has  been  determined,  the 
Strehl  relation  is  used  to  compute  the  intensity  degradation 
due  to  thermal  blooming. 

i  =  exp  (-a2)  (2.  13) 

T  o 

The  result  of  equation  2.13  is  a  relative  intensity  (Irel) 
ratio.  To  is  the  ideal  on-axis  irradiance  with  no  phase 
distortion.  The  Strehl  relation  above  is  thought  to  be  too 
severe  a  model  for  Irel  below  0.3  [Ref.  6].  For  that 
reason,  if  a2  is  less  than  1.2  ,  the  Irel  will  oe  computed 
using  polynomial  curve  fits  developed  from  GUTSMTF  results. 
GUTSMTF  is  a  full  wave  optics  code  utilizing  fast  Fourier 
transforms.  For  a  description  of  the  curve  fit  method 
above,  see  the  subroutine  BLOOM  explanation.  [Ref.  7] 

Combining  the  thermal  blooming  effect  with  the  other 
effects,  such  as  turbulence  and  jitter,  is  accomplished  by 
averaging  the  results  of  two  different  approaches.  The  first 
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approach  is  the  ESS  (root  sum  squared)  method.  This  method 
cx  combining  tne  Ire 1  due  to  the  effects  of  thermal  ticcming 
with  the  Irel  due  to  gitter  ,  team  quality,  and  turhulence  is 
accomplished  as  follcks 


Irel 


rss 


(2.  14) 


where 
res  ul 
plica 


The  t 
aver  a 
atten 


tases 
and  i 
methc 
te  to 
misti 
of  tai 
full 

1 

used 

menti 


irel  ,  is  the  thermal  ticcming  result  and  Irel  is  the 
tb  a  o 

t  of  the  other  effects.  Ihe  second  approach  is  a  multi- 
tive  approach  arc  is  simply 


Irelm  "  Irelo  X  Ireltb  <2-15> 


kc  ccmtined  Irels  cttained  ty  these  methods  are  then 
ced  tc  give  the  total  intensity  ratio  due  to  all  the 
uating  or  distcrtiny  propagation  effects. 


I 

Io 


Irel 


tot 


Irel  +  irel 

m  rss 

2 


(2.  16) 


the  ideal  diffraction  limited  on-aris  intensity.  Ihe 
for  the  above  averaging  process  is  empirical  in  nature 
s  ar  attempt  to  adjust  the  results  obtained  ty  the  E IS 
c  alone.  The  results  produced  by  HSS  were  thought  tc 
c  optimistic.  Ihe  multiplicative  method,  a  mere  pessi- 
c  approach,  was  therefore  included.  Ihe  ultimate  Ireis 
ted  are  very  clcse  to  these  obtained  ty  GUTSMTF,  the 
lave  cede.  [Bef.8] 

here  are  some  limitations  tc  the  thermal  blooming  model 
in  GC1SAVG  in  addition  tc  the  assumptions  already 
cned.  First,  the  wind  (V0)  is  applied  as  a  corstant 
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C.  SCATTERING  AND  APSCRPTICN 

GETSAT7G  does  not  contain  ariy  type  ox  atmospheric  mccel 
with  respect  to  scattering  and  absorption  data.  Extinction 
coel ficients  must  be  entered  by  the  user  in  the  appropriate 
subroutines,  ALPS  and  ALPA,  or  the  program  can  be  modified 
to  accept  a  separate  cata  file.  Once  Inserted  into  the 
program,  scattering  ard  absorption  coefficient's  are  treated 
without  distinction  between  aerosol  and  molecular  mecha¬ 
nisms.  Therefore,  the  coefficients  used  must  represent  the 
total  effect  of  scattering  or  absorption. 

a(h)  =  a(h)  .  +  a(h)  i2m  17) 

mol  aer 


ie 


o(h) 


(2-  18) 


Tran  smi  ssi  on  due  tc  scattering  and  absorption  are 
computed  indentically  and  given  by 


•(2.  IS) 


where  0  is  the  zenith  angle,  ard  the  integration  limits  are 
the  altitude  of  the  transmitter  and  the  vertical  extent  or 
the  atmosphere.  The  extent  of  the  atmosphere  in  the  program 
is  defined  as  30  km.  Figure  2.3  snows  the  program  applica¬ 
tion  cf  scattering  ard  absorption. 

E.  E  EAH  COflLITI 

GllSATiG  allows  the  user  tc  describe  the  neam  at  the 
aperture  in  terms  of  the  electromagnetic  field  amplitude 
distribution  and  total  beam  yuality.  The  intitai  amplitude 
field  used  by  the  program  is  Gaussian  in  shape  with  the  user 
specifying  e  waist  diameter.  The  waist  diameter  is  defined 
by  the  1/e2  point  on  the  distribution.  Truncation  cf  the 
this  Gaussian  field  will  of  course  depend  on  the  aperture 
diameter  and  the  size  of  the  central  obscuration.  If  the 
waist  diameter  is  made  large  compared  to  the  aperture,  a 
more  uniform  distribution  results. 

Eeam  guality  at  the  aperture  exit  may  be  specified  by 
two  different  parameters.  One  cf  these  is  the  'times 
diffraction  limited  number' ,  N.  N  has  been  used  in  a 
general  way  to  mean  at  increase  in  far-field  spot  size  or  as 
a  ' pcwer-in-the-bucket'  ratio.  In  GUTAVG,  N  is  a  total  team 
guality  term.  The  second  beam  guality  parameter  is  a  ncndi- 
uensicnal  term  representing  the  BMS  phase  distortion  at  tie 
laser  wavelength  at  the  exit  aperture,  ^rms  - 


X 


Ihe  parameters  axe  related  to  each  ctner  and  to  the  inten¬ 
sity  degradation  by 


I 

Io 


1 

N2 


/  2  7T  6  \ 

exp  -  rras  ) 

V  a  / 


(2.20) 


1c  apply  beam  quality  tc  the  propagation  problem,  ar  hlF 
array  is  constructed  representing  a  phase  screen  at  the 
aperture.  He  phase  is  assumed  tc  oe  a  Gaussian  random  vari¬ 
able  kith  a  zero  mean  value.-  For  the  axi-symmetric  team 
coasicexed  in  GU1AVG,  the  dlF  is  [Ref.  11] 


Mb(p) 


=  exp  ^-k  2|cr 2  —  C ^  (  p ) J 


(2.21) 


CA  is  the  autocorrelation  function  of  the  phase  and  is 


aerrnea  as 


C^(p)  =  a  2  exp 


-(f)' 


(2. 22) 


2tt6 


where 


is  the  phase  variance,  (■ 


rms  2  . 

— -  )  ,  and  L  is  the 

phase  ccrxelation  length  £Ref.  12].  The  neam  quality  a  IF 
array  is  combined  with  th,e  MTF  arrays  due  to  ether  propaga¬ 
tion  effects  to  determine  the  complete  system  MTF  and, 
hence,  tie  total  irradiance  degradation.  Based  on  the 
provided  input  parameter.  Figure  2.4  shows  the  general 
treatment  cf  beam  guality  within  the  program. 


1CBECIINCE 

GtlSAVG  uses  the  Eufnagel  model  £Re£.  13]  for  as  an 


indicator  of  the  optical  effects  of  turbulence  along  the 
propagation  path.  is  the  refractive  index  structure 
constant  and  represents  the  refractive  index  in  the 
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atmosphere  as  a  function  of  turbulence  induced  density 

fluctuations.  Ihe  model  is  an  empirically  derivea  verc 

profile  cf  C2  . 

n 

Iriec  [Bef.  14]  has  developed  a  parameter  which  is 
directly  related  to  tie  tehavicr  cf  a  coherent  bean  in 
turbulent  medium.  This  term  is  called  the  effective  cc 

ence  diameter,  r  .  Ir  the  case  of  a  laser  transmiter, 

o 

represents  a  physical  limit  to  the  transmitter  diamete 
near  diffraction-limited  beam,  for  a  transmitter  diaae 


larcer  than  r 


degradation  cf  the  ream  by  turtuler 


will  cccir.  If  D  is  smaller  than  r, 


tnen  near  diifra 


limited  propagation  will  be  achieved.  Beam  wander  or  p 
tilt  occurs  for  transmitter  diameters  approximately  eg 
rQ.  lura  [Befo  t5  ]  has  defined  a  somewnat  different  L 
related  term  that  car  be  thought  of  as  a  lateral  coher 
length,  p0.  These  twc  guantities  are  given  by 


r  u  = 


fri!  k2sec0 
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where  k  =  \  ,  0  is  the  zenith  angle,  and  the  limits 

integration  are  b  ,  the  vertical  extent  of  the  atmcs 

atm 

and  tt  ,the  height  of  the  transmitter. 
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Po  may  be  specified  by  the  user  as  a  program  input  or 
the  program  may  be  allowed  to  compute  it  based  oi  equation 
2.2U  .  The  user  may  also  input  Po  indirectly  by  specifying 
the  ’seeing  conditions’  ,  a  quantity  used  by  astronomers  to 
describe  angular  spread  of  a  stellar  point  source  (see  input 
section  and  Figure  2.5).  I  he  program  uses  P0  to  compute  the 
atmospheric  MIF. 

The  HTF  of  the  turbulent  atmosphere  is  deteris  ined  by 
developing  the  structure  function  of  the  turbulence.  This 
development  is  demonstrated  by  Yura  [Ref.  16].  The 
resulting  MTF  is  given  by 


(2.  26) 


This  HTF  effectively  applies  the  the  effect  of  tarbulence 
along  the  propagation  path  as  a  phase  screen  at  the 
a pert  ure . 

F.  JITTER 

Beam  jitter  is  a  user  input  and  is  specified  by  the 
variance  of  the  angular  excursions  of  the  beam.  J  sing  the 
2-siyma-p  definition. 


(2.27) 


=  V2  ( a2  +a2 ) 
x  y 


2a 


P 


where  px  and  Py  are  the  axial  variances  of  the  jittered  beam 
center  in  the  far-field.  px  and  Py  are  random  variables  with 
Gaussian  distribution  and  in  the  symmetric  case  ,  as  consid¬ 
ered  by  GUTSAVG,  px=  Py  .  The  resultant  intensity  distribu¬ 
tion  due  to  jitter  will  also  be  a  Gaussian  distribution  with 
(20p)  representing  the  spot  radius  defined  at  the  1/e2 
point.  In  other  words,  86.5%  of  the  beam  energy  will  reside 
within  the  radius  2a  .  [Ref.  17] 
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Jitter  can  be  skvn  to  he  a  wavefront  tilt  at  the  aper¬ 
ture.  Using  this  approach,  a  phase  screen  at  the  aperture 
can  he  used  to  characterize  the  effects  of  jitter  and  a 
jitter  Ml?  developed.  That  MTF  is  given  by 


ML  (  p)  =  exp 


-k2  o2  (2op)  2 
8 


(2.28) 


figure  2.6  is  a  general  flow  diagram  for  the  treatment  cf 
jitter  within  tiie  pregram. 


G.  AIAE1IVE  OPTICS 

The  ability  to  apply  adaptive  optics  corrections  to  the 
propagation  problem  has  bee n 'included  in  GUTSAVG.  The  user 
has  several  options  kith  respect  to  the  type  and  degree  or 
compensation  desired.  The  following  general  discussion  and 
figure  2.8  provides  the  needed  insight  to  the  effects  of 
selecting  the  adaptive  optics  options. 

The  najcr  compensation  mode  provided  by  the  program  is 
invoked  by  selecting  full  zonal  adaptive  optics  with  consid¬ 
eration  given  to  isoplanatic  effects.  When  selected,  this 
model  results  in  the  correction  of  beam  degradation  due  to 
turbulence.  This  is  accomplished  by  correcting  turbulence 
induced  tilt  and  then  adjusting  Po  so  as  to  produce  a 
predetermined  level  cf  adaptive  optics  performance.  This 
predetermined  performance  is  as  measured  by  the  Strehl  ratio 
given  a  residual  phase  variance  determined  by  the  adaptive 
optics  sensor  phase.  Parameters  determining  the  phase  error 
are  the  response  bandwidth  of  the  adaptive  optics  system, 
the  number  cf  system  actuaters,  the  reflected  radiant  inten¬ 
sity  of  the  target,  and  the  target-to-sensor  transmission . 
The  resultant  .Po  found  in  this  manner  is  then  used  to 
compute  the  atmospheric  MTF.  When  this  ’corrected’  MTF  is 
used  to  compute  the  far-field  intensity,  the  result  will 
represent  an  adaptive  cptics  corrected  value. 
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Without  invoking  full  zciia  1  adaptive  optics,  the  user 
cay  apply  a  tilt-only  correction  for  turbulence.  In  this 
case  some  or  all  of  tie  tilt  due  to  turrulence  may  be 
removed  before  computing  the  atmospnenc  rilr.  The  degree  cf 
tilt  compensation  is  specified  by  the  user. 

As  mentioned  above,  isoplaratic  effects  are  included  in 
the  adaptive  optics  calcula tic  ns .  This  is  also  an  option, 
however,  and  isoplanatic  calculations  may  be  inhibited  by 
the  user.  Tbe  effect  considered  is  the  limitation  cf  the 
adaptive  optics  system  given  an  isoplanatic  angle  smaller 
than  the  target  lead  angle.  Fried  pfief.  18]  provides  a 
discussion  cf  iscpiaratism  and  development  cf  the  iscpla- 
ratic  arcle. 

Although  adaptive  optics  compensation  for  thermal 
blooming  is  not  modeled  in  a  strict  sense,  a  thermal 
blooming  correction  factor  can  be  applied.  'Inis  factor  is 
simply  a  fractional  constant  that  multiplies  the  thermal 
blooming  phase  variance  before  the  Strehl  relation  is  used 
to  compute  the  intensity  degradation. 

Eefer  tc  the  input  definitions  and  the  subroutines 
involved  with  adaptive  optics.  In  addition,  refer  tc  figure 
2.7  fci  more  explanation. 
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Figure  2.1  GC1SAVG  Ihermal  Blooming  Algorithm. 
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Figure  2«2  Tiermal  Elccming  Algorithm 
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(coat) 


Figure  2-3  GU2SAVG  Scatteriug  and  Absorption  Algorithm. 
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Figure  2.4 


G01SAVG 


Eeam  Quality  Algorithm. 
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figure  2.5  GOISAVG  lurfculence  Algorithm. 
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figure  2.i 


GJJTSAVG  Jitter  Algorithm. 
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INPUT  DECISIONS 


Nf-LAGA  FULL  ZONAL  AO 
1  YES 

0  NO 

NO  I  SO  INHIBIT  ISOPLANATIC  CALCULATIONS 

I  YES 

0  NO 

AOAP  RESIDUAL  TURBULENCE  TUT 

0.0  REMOVE  1001  OT  THE  TURBULENCE  TILT 
1.0  REMOVE  0%  Of  THE  TURBULENCE  TILT 


figur€  2.7  GZIS&IG  Adaptive  optics  Algorithm. 
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F i  gore  2.8  GOTSA\G  Adaptive  Optics  Algorithm  (ccot). 


Figure  2.9  GUTS  AUG  Adaptive  Optics  Algorithm  (cent) 


Figure  2.10 


GOTSADG  Adaptive  Optics  Algorithm  (cent). 
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Fig  ure  2.11 


GUTSAiG  Adaptive  Optics 


Algorith* 


(cent) . 


35 


Figure  2..  12 


GUTS AVG  Adaptive 


Optics  Algorithm 


(cent) „ 
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III.  PBOGBAJ  AND  SDEP BQGBAfl  DZSCBIP1ION 


A.  I BCGIAfi  INPUTS 

lie  following  are  tie  inputs  necessary  to  utilize 
GUTSATiG.  Default  values  are  indicated  for  parameters  not 
absolutely  required  for  program  operation.  At  NPS  these 
input  parameters  are  entered  via  an  input  file.  A  copy  or 
this  file  is  provided  in  Appendix  A.  (*)  indicates  a  ncndi- 
mensicral  parameter. 

•  DIA  (meters) 

lie  diameter  of  tie  transmitter  aperture. 

•  IIACE  £  (meters) 

lie  diameter  of  the  central  obscuration  of  the  trans¬ 
mitter. 

•  EIAMS2  (meters) 

lie  Gaussian  waist  diameter  of  the  amplitude  distribu¬ 
tion  at  the  aperture,  (measured  at  the  1/e2  point) 

«  i A Tii  (meters) 

Ihe  laser  wavelength. 

•  PICIAL  (Watts) 

lie  total  power  at  the  aperture. 
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«  1EIIH1  (*) 

Ihis  is  the  cftei  used  "tines  diffraction  limited 
number"  and  represents  total  team  quality.  As  used  in 
GC1S,  it  is  related  to  the  BUS  phase  distortion  at  the 
aperture  by  the  Strehl  approximation 


(TDFLMT) 


2=  exp 


2  it  6  '  z 

rms 

A 


(3.1) 


(1EF1M1)2  ,  therefore,  is  equivalent  to  the  ratio  of  the 
cr-axis  diffraction  limited  intensity  to  the  cn-axis 
intensity  resulting  from  the  near-field  phase  distor¬ 
tion,  6 


rms 


(TDFLMT) 


I 

I  o 


(3  .2) 


For  a  discussion  cr  the  limitation  of  3.1,  see 
£  Bef  .  IS]- 

•  RAVEEC  (*) 

Ihis  term  is  the  E MS  phase  distortion  at  the  laser  aper¬ 
ture,  ncndi mensicnalized  hy  the  wavelength. 


WAVEBQ  = 


rms 


(3.3) 


•  SCA1EQ  (meters)  default  =  DIA/o 

Ihis  term  is  the  transverse  phase  correlation  length  at 
the  aperture. 
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THSEE  (arcsec) 


default  =  f(RHOO) 


This  is  a  qualitative  term  used  by  astronomers  to 
describe  'seeing  conditions'  in  the  visible  range.  If  a 
point  source  is  viewed  from  the  earth,  it  may  not  appear 
as  a  point  source  but  as  a  'smear'  or  spot.  The  angular 
spread  of  this  spot  is  the  parameter  THSEE.  If  RHOO  is 
not  specified  as  an  input  and  THSEE  is,  the  program  will 
use  THSEE  to  compute  RHOO.  Other  than  this  case,  THSEE 
is  not  used  but  will  be  computed  as  a  function  of  RHOO 
and  included  as  an  output. 


•  HGRND  (meters) 

The  height  of  the  ground  at  tne  transmitter  position 
above  MSL. 

•  H TRANS  (meters) 

The  height  of  the  transmitter  above. 


•  HS  AT  (aeters) 

The  orbital  altitude  of  the  satellite  above  1  SL  (at 
zenith)  . 


•  THETHX  (degrees) 

THETMX  is  the  angle  measured  from  zenith  below  which  the 
laser  will  not  transmit  for  a  zenith  pass.  For  offset 
flight  paths,  it  should  be  noted  that  the  transmitter 
will  point  below  this  value.  For  an  illustration  of  the 
engagement  envelope,  see  figure  3.1 
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ICEF  (meters) 


Per  a  target  that  dees  ret  pass  direct  overhead,  this 
input  specifies  tie  amount  that  the  target  ground  track 
is  offset  from  the  ground  track  of  the  overhead  case.  It 
is  the  distance,  as  measured  from  the  transmitter,  to  a 
per p endiculaz  intersection  cn  the  ground  track. 


■»  EHCO  (meters)  default  =  f(IriSEr)  or  aHCISE 

Ibis  is  the  turbulence  coherence  length  as  defined  by  H. 
I-  Yura.  £Eef.  2C  ]-  (see  turbulence  section  for  mere 
discussion  ard  references) 


•  VO  (m/sec) 

He  atmospheric  wind.  The  direction  of  the  wind  is 
parallel  and  opposite  to  the  direction  of  taryet  metier, 
licte,  in  the  present  program,  VO  is  a  constant  indepen¬ 
dent  of  altitude. 

•  ACIICfl  (*)  0.0  to  1.0 

This  parameter  allows  tie  user  to  correct  for  thermal 
ilccming  as  if  by  adaptive  optics.  The  value  entered  may 
range  from  0.0  tc  1.0.  If  0.0  is  used, complet e  thermal 
tlccairg  compensation  will  occur.  Conversely,  if  1.0  is 
entered,  no  compensation  will  be  applied.  The  variance 
of  the  phase  distortion  is  multiplied  by  this  correction 
factor  ref ore  the  Strehl  relation  is  used  to  compute  the 
relative  intensity  reduction  due  to  thermal  blooming. 

Y  =  exp ( -O2  *  AOBLOM)  (3.4) 

o 
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•  MGS PI  (*)  0..C  to  1.0 


AVGSPT  allows  th€  user  latitude  in  defining  the  :ar- 
field  spot  size  tc  otner  than  that  indicated  by  GUTSAVG 
analysis.  A  value  of  1.0  would  result  in  the  peak  irra- 
diance  according  to  the  program  analysis.  An  entered 
value  of  .5/  for  example,  would  result  in  a  peak  lrradi- 
ance  50S  less  than  the  program  would  otherwise  indicate. 
AVGSET,  then,  is  an  adjustment  factor  that  allows  the 
user  tc  account  for  effects  not  addressed  in  the  G01SAVG 
propagation  calc ulations. 

•  SIGJI1  (radians) 

SIGJI1  is  the  2o^  variance  for  pointing  and  tracking 
jitter. 


*  AC  AP  (*)  0.  C  to  1.0 


Ihis  te 

rm  is  a  correction 

factor 

for  tilt  due 

to  turbu- 

lence. 

It  represents  the 

residua 

1  tilt 

after 

AO  compen- 

saticn. 

ADAP  may  be  varie 

d  from 

0.0  to 

• 

O 

1 

value  cf 

1.0  would  result  in  no  tilt  due  to  turbulence  being 
removed  while  a  value  of  0.0  would  result  in  total  tilt 
due  to  turbulence  'compe nsation.  If  full  AO  is  selected 
the  program  will  set  ADAP  equal  to  0.0. 

•  NF1AGA  (*}  0  cr  1 

NFIAGA  is  a  selection  indicator  for  full  zonal  adaptive 
optics.  Enter  a  0  if  AC  is  not  desired,  enter  a  1  if  AO 
is  desired.  If  AC  is  used,  XJT,BWIDTH,  and  NA  must  be 
specified. 


4  1 


O  HC3SC 


(*)  0  cr  1 


HCISC  is  a  selection  indicator  for  isoplanatic  calcula¬ 
tions.  Enter  a  1  to  inhibit  isoplanatic  caiculaticrs; 
enter  a  0  if  isoplanatic. calculations  are  desired. 

•  XJ1  (Watts/sterad) 

XJ1  is  the  target  radiant  intensity.  It  is  one  of  the 
factors  used  to  determine  how  noise  arfects  the  response 
cf  the  AO  system. 

«  ElilElE  (Hz) 

EHIE3H  is  the  bandwidth  of  the  adaptive  optics  system- 

•  HA  (*) 

HA  is  the  number  cf  AO  actuaters  used  tc  perform  phase 
ad jrstment . 

•  AEE1C2  (*) 

This  is  the  percent  tra nsmissioh  at  zenith  at  the 
sensing  wavelength  of  the  AO  system.  This  term  is  used 
in  the  determination  of  relative  noise  at  the  AO  sensor. 

•  HI  (*) 

The  rusher  of  iteration  steps  for  the  angle  loop  cf  the 
program  from  THEIKX  to  0. 
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•  N 2  (*) 

lie  number  of  altitude  intervals  for  absorption  and 
scattering  determination. 

«  S3  (*) 

lie  number  of  altitude  intervals  for  the  turbulence 
calculations. 


*  N  4  {*) 

lie  number  of  iteration  intervals  rcr  the  .1TF  calcula- 
tic  ns. 

•  S5  (*) 

lie  number  of  intervals  for  the  slant  path  update  c± 
thermal  blooming. 


£.  ENGAGEMENT  GEOHElfT 

lhe  target  engagement  -geometry  is  that  of  a  earth-tased 
transmitter  and  a  target  satellite  at  a  given  orbit.  Nc 
attempt  is  made  to  define  an  engagement  envelope  based  cn 
threshold  irradiance  cr  fluence1.  The  input  parameters 
defining  the  engagement  window  in  GUTSAVG  are  THET2JX,  1CFF, 
and  HSAI.  lie  general  geometry  of  this  window  is  shewn  in 
Figure  3.1  .  Only  half  of  the  total  transit  window  is 
addressed  in  the  program  calculations ;  the  results  are  the 


1  lhe  GUISFP  (footprint)  version  of  gets  was  writter 
to  dc  this.  Except  for  this  feature,  the  propagation  calcu¬ 
lations  are  identical  to  GUTSAVG. 
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came  fcr  either  naif.  The  program  output  reflects  this 
half-window  evaluaticr  except  for  parameters  such  as  total 
fluence  and  shot  time  which  are  simply  double  the  computed 
values. 

figure  1.2  defires  the  earth  center  ancle  (BCANG) .  ECANG 
is  a  function  of  the  user  input  1H2TMX.  Most  of  the  geome¬ 
tric  calculations  are  referenced  to  earth  center.  Therefore/ 
this  argle  is  used  fcr  computing  such  positioning  data  as 
the  angle  interval  at  which  the  irradiance  will  oe  evaluated 
(see  f  i  gure  3.3)  . 

Offset  flight  paths  reguire  a  coordinate  translation  as 
snown  in  Figure  3.4  .  Position  and  velocity  relative  to  the 
transmitter  are  computed  as  in  Figure  3.5  .  It  should  he 
noted  that  if  the  flight  path  is  offset,  the  zenith  angle 
will  exceed  IHETMX  fcr  part  or  all  of  the  window.  This  is 
because  tie  window  is  defined  in  the  x-z  plane  only. 


44 


Zenith 

Flight 


Hs  -  Zenith  altitude  of  target 

0^  -  Maximum  zenith  angle  for  engagement  (zenith  path) 

L  -  Offset  distance 
o 


Figure  3-1  General  Engagement  Geometry. 
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z 


figure  3.2  Earth  Center  Angle  (ECANG)  . 
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I  CM 


z 


Figure  3.3  Angle  Intervals  for  Irradiance  Evaluation. 
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s 


figure  3.4  Coordinate  Translation  due  to  Plight  Path  Cfiset. 


4  £ 


Velocity  components 


z 


V 

x 

V 

y 

v 

z 


=  -V  cos  (0  ) 
s  e 

=  V  s in ( 6  ) cos  (y) 
s  e 

=  V  sin ( 6  ) sin (y) 
s  e 


V 


Position 

X  =  R  sin ( 0  ) 
s  e 

Y  =  Ly-cos (y) + (Z 0-R-Lz) 

• s in (y ) 

Z  =  -Ly ■ sin ( y) + ( Z o -R-Lz ) 
•  cos ( y) 

where  Z0=  R  cos  (0  ) 
u  s  e 


Range  from  transmitter (T)  to 

satellite (S ) 

R,  =  Z  +  X  +  Y 
ts 


(For  Ly  and  Lz  definition,  see  Figure  3.4) 


figure  3.5  Velocity  and  Position  Relative  to  Transni t ter. 

4  S 


C.  BAIN  EBCGBAB  FlOfi  DIAGBAB 

lie  mai r  program  ircstiy  consists  of  geometry  calcula¬ 
tions  and  decision  flew  points.  The  decision  points  allcw 
branching  tc  adaptive  optics  aDd  isopianatic  subroutines. 
Bost  cf  tie  propagation  calculations  are  done  within  the 
subroutines.  Only  tve  major  iteration  loops  reside  in  the 
main  program,  the  ancle  interval  loop  and  the  combined  Klf 
integraticr  loop.  Tie  following  is  a  general  flow  diagram 
for  tie  nair  program.  Decision  variables  are  user  inputs  cr 
program  defined  parameters. 
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figure  3.6  iiain  Pregram  Flow  Diagram. 

5  1 


Figure  3.7 


Bain  Program  Flow  Diagram  (coot) 


Bain  Program  Flow  Diagram  (cont) 


Figure  3.8 


Figure  3„ 9 


Bain  Program  Flow  Diagram  (cont) 


Figure  3-10  Bain  Pregram  Flow  Diagram 


=;  c 


(cont) 


eain  Program  Flow  Diagram  (contj . 
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Figure  3U 1 1 


Figure  3- 1 2 


Bain  Program  Flow  Diagram 


(cont)  . 
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0 


compute:  isoplanauc 
angle  ANO  test,  based 
on  look-aheao  angle, 
to  determine  if 

FULL  ISOPLANATIC 
CALCULATIONS  need  to 
BE  DONE 


NISO-1  YES 
NISO-O  NO 


COMPUTE  IPRAOIANCE 
REDUCTION  OUE  TO  SLANT 
RANGE. 


TLTISO 


COMPUTE  WAVEFRONT  TILT 
OUE  TO  ISOPLANATIC 
EFFECT 


figure  3o 1 3 


lain  Program  flow  Diagram  (cont) • 


se 


figure  3-14  Hair  Program  flow  Diagram  (coiit). 


5  S 


Figure  3.15 


ilaia  Program  Flow  Diagram  (cont)  . 


6  C 


figure  3. 1 6 


Bairn  Program  Flow  Diagram 


(CODt)  . 


6  1 


L 


£  CE-ECOHNS  DISCfi If  T1CN 


1.  AE5CR3  aid  5 CAT 


As  noted  in  the  preceding  chapter,  absorp  tier  and 
scattering  effects  axe  treated  indent ically .  The  AESCRE  axd 
SCAT  algorithms  perform  the  integration 


T 

a 


T 

s 


exp 


exp 


h 


atm 
a (h) dh 


atm 
a (h) dh 


(3.5) 


(3.o) 


The  result  is  the  transmission  due  to  the  rotai  absorption 
cr  total  scattering.  The  correction  for  slant  range  is 
applied  kithin  the  angle  interval  loop  ox  the  main  program 
and  is  siaply 

( sec  0 ) 

T  =  T  (3.7) 


is  the  xerith  angle  of  the  target  with  respect  tc  the 
transmitter . 

The  extinction  coefficients  for  absorption  (a)  arc 
for  scattering  (o)  axe  provided  to  the  subroutines  by  a  call 
to  routines  ALFA  and  A.1FS  respectively. 
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r 


TABU  I 

AESCfiB  and  SCAT  Program  Variable  Definitions 


Vjriable 

Fortran 

III - 

name 
~ 3_ain 

Definition 

1  (AISCEB) 

ci 

T 

TABSC 

Total  molecu lar 
transmission 

Is  (SCAI) 

T 

TSCAIC 

Total  scattered 
transmission 

h  (neters) 
atm 

HATHC 

HATMC 

Height  of 
atmosphere 

ht  (meters) 

HT 

H TEA  S3 

Height  above  MSI 
of  transmitter 

a  (h)  (km)  1 

ALP 

- 

Ansorpticn  coef. 
at  given  altitude 

a  (t)  <kn)  1 

A1S, 

- 

Scattering  coef 
at  given  altitude 

- 

H 

N  2 

#  of  integration 
intervals 

2-  A1FS  and  AIFA 


These 

two  routines 

are 

also  i 

dentical. 

Th  ey 

p r c vide 

the  extinction 

coefficients 

tc 

ABSOEB 

and  SCAT 

for 

a  speci- 

fied  altitude. 

At  each  alti 

tude 

where 

a  coefficient 

value  is 

desired,  a  linear  interpolation  is  performed  between  data 
points  supplied  by  tie  user.  Data  statements  precede  each  of 
these  subprograms,  and  it  is  with  these  statements  that  the 
absorption  and  scattering  data  should  be  entered.  Units  for 
the  coefficient  and  corresponding  altitude  should  be  km  1  and 
km  respectively. 
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T  ABIE  II 


AXIS  and  A IF  A  Program  Variables  and  Definitions 


Variable 

Fortran 

511 - 

name 

~Hann 

Definition 

a(t) 

(kc)  1 

S 

A  IS 

Scattering  ooef. 
at  specified 
altitude 

a(k] 

(bn)”1 

A 

Aif 

Absorption  coef. 
at  specified 
altitude 

h  (meters) 

H 

- 

specific  altitude 

— 

AIT  (FI) 

- 

Altitude  data 
list 

— 

ATA  (SI) 

- 

Absorption  data 
list 

— 

ATS  (M) 

- 

Scattering  data 
Inst 

NI 

dumber  of  points 
in  data  list 

COCST 

Fen  a  gaussiar. 
amplitude  distribution 


beam  with  constant 
at  the  transmitter 


phase,  the  initial 
apertuc  e  is 


U  (r) 


U  0exp 


(r/w) 


(3.8) 


where 

UQis 

the 

ampli 

ttde  and  w 

The  p 

ur  pose 

of 

HOC  ST 

is  to  com 

given 

aper  t 

ure 

pewe  r 

m 

Ey 

the 

sea  la 

r  wave  app 

distr 

itutic 

r  is 

:  [Ref 

.  22] 

is  the  spot  size  I  Ref. 
pute  the  constant  U0  for 

nomination,  the  inters! 


t 


21]. 

a 


y 


64 


I  (r)  =  | U ( r )  |  2 


(3.9) 


lo  relate  tie  field  distribution  to  the  power,  toe  intensity 
distribution  is  intecrated  over  the  aperture. 


,t  -  ,  f 


I (r) dr  2 


(3w  10) 


Substituting  3.8  and  3-9  into  8-10 


-P 

IT  . 

1 


Pt  =v  I  U0  exp 


- (r/w) 


dr  : 


(3.  11) 


then  rearranging  terns,  produces  tne  expression  for  U  . 


<3.  12) 


Ihe  intecratior  limits  r  and  r.'  are  the  radius  of  the 

o  1 

transnitter  and  the  obscuration,  respectively. 

LOCST  perforns  the  integration  in  3.12  using  the 
trapezoidal  rule.  A  call  is  made  to  subroutine  FIELD  to 
evaluate  exp[-(r/w)2]  which  is  simply  U  (r)  with  U0=  1.  For 
this  reason,  U0is  defined  as  unity  when  FIELD  is  called  by 
U0CS1 . 
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TA  f IE  Ill  J 


DOCS  I 

Pr  o  g  r  a  a 

Variables 

and  Definitions 

Variable 

For  tran 
3uE - 

lid  E6 

"lain 

Definition 

*5  “  1 

B  o  (*)  E 

uo 

- 

Normaiizatio  n 
constant 

It  (Katts) 

p 

P1C1M 

Total  aperture 
exit  power 

- 

N 

ii  1 

n  of  inteora tioc 
increments 

rQiieteis) 

SO 

- 

Outer  radius  cf 
transmitter 

r  j_  {mete  rs) 

F.I 

l 

UB 

- 

Radius  of 
obscuration 

0  (r  )  i'n  )  2m 

- 

Field  ampiit ude 
at  radius  r 

v  (meters) 

BUB  A I 

- 

Radius  of  spot 
at  aperture 

a  *  III ID 

f  HID  calculates  tie  field  distribution  f  or  a  axi- 
symmetric  Gaussian  beam  with  constant  phase. 


U (r )  =  U  o exp 


- (r/w) 


(3. 13) 
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TABU  IV  j 

FIELD  Program  Variable  Definitions  1 


1 

4 

i 

[ 

Variable 

h 

Fortran  name 

I’ll  ~TI  ain 

Definition 

« 

i 

1 

i 

Oflr)  ( K )  mn 

k 

UR 

Field  amplit 
at  r 

ud€ 

i 

l 

i 

1 

s» 

o 

IZ) 

UO 

Nor malizatio 
constant 

n 

1 

i 

i 

w  (meters) 

3MEAI 

Radius  of  sp 
at  aperture 

ot 

- •  LAB  PIP 

In  the  Fraunhofer  region,  the  amplitude  distribution 
can  be  found  by  taking  the  Fourier  transform  of  the  aperture 
distribution  [Ref.  23].  Using  polar  coordinates  and  noting 
axial  symmetry,  U(p)  in  the  far-field  can  be  expressed  as 


2  7T 


U(p)=  Xi-  / /u<r>expr-(i?2) 

a  0  L  \  / 


(i  2tt. 

-y-H  •  rp  •  cos0 


rdrdG 


(3.  14) 


tsing  integral  relation  for  the  zero  order  Bessel  function, 
equation  3.14  can  be  simplified  to  [Ref.  24] 


U(p)  = 


2tt 

Az 


/  J.(^) 


rdr 


(3.  15) 
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lie  intensity  in  the  far-field  is  then  given  by  I  (r)  =  |U{r)|2 
cr  in  tens  cf  eguaticn  3.15 


(3.  16) 


where  i  and  a  are  the  radius  ci  the  transmitter  and  the 
radius  ci  the  obscuration,  respectively.  [Ref.  25] 


IJSRilD  evaluates  3.16  at  a  specified  number  cf 


increments  in  the  fan-field  and  assigns  these  values  tc  an 
array  F  (i)  .  F{i)  is  normalized  with  the  total  aperture  power 
so  that  the  HTF  prod  coed  would  be  unity  at  the  origin,  hcte, 
however,  that  when  tie  MTF  is  computed  in  this  case,  the  iv.IF 
is  normalized  so  as  tc  produce  an  Irel  {intensity  relative) 
value  when  integrated. 


r 


0 


I  0  is  determined  by  evaluating  eguation  3.16  with  r  =  C.O. 
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TABLE  V 

, 

f  ABFLD 

Program 

Variables 

and  Definitions 

Variable 

Fortran 
Suf - 

name 
~H  am 

Definition 

k  (m)"1 

CKR 

- 

(  2tt/A) 

P 

El 

- 

fiadius  in 
f ar-f ield 

r 

EO 

- 

Radius  in 
aperture 

I0  (Katts) 

PMAXC 

PMAXG 

Far-fiela  oQ-axis 
intensity  for  a 
dif fr action 
limited  team 

Id)  (Vi/m2) 

F  (I) 

TISC 

Irradiance  in  , 

far-f ield 

Ha(p) 

G  (I) 

IEEMF 

Aperture  MTF 

- 

DX2 

DX 

Aperture  MTF 
increments  , 

E  t(Hatts) 

P 

P1C1AL 

Total  power  in 
the  aperture 

_ 1 

e .  EGIFEL 

If  the  'times  diffraction  limited  number*  (TDFLMT)  is 
not  specified,  this  routine  computes  the  Irel  value  due  tc 
team  guality  which  is  then  used  to  compute  TDFLMT  in  the 
main  program.  The  calculation  performed  is 


Irel 


2tt 


Ma(p)  Mj^(p)  pdp 


(3.  18) 


where  M 
guali ty 
that  th 


ais  the  aperture  MTF 
phase  screen.  Ma  i 
is  integration  produ 


and  M^is  the  MTF  for 
s  defined  in  subroutin 
ces  an  Irel  value.  M £ 


the  beam 
e  FAEFID  so 
is  the  team 
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guality  K1F  as  computed  by  subroutine  alFBQ.  1DFL Ml  is 
commuted  in  the  main  program  aDd  is  1/  Irel  . 


TABIE  VI 

ECIHEL 

Pro  cram 

Variables 

and  Definitions 

Variable 

Fortran 
Sul - 

na  me 
~HaIn 

Definition 

K 

— 

TEF1B1 

limes  diff. 
limited  # 

Irel 

EEL 

TEQ 

- 

ea(M 

A(I) 

IEEMF 

Aperture  MTF 

»b<e> 

F{I) 

— 

MTF  for  phase 
screen 

- 

N4 

N  4 

#  of  iterations 
for  iTIF  calcula  tic  ns 

I 

l. 


1 .  IEJUAV 


If  the 
rot  specified, 
approximation 


E  MS  phase  distortion  parameter  (WAVEEQ 
this  routine  computes  it  based  on  the 


is 


Irel 


exp  (-a2) 


<3.  19) 


2  it  6 

where  a  =  — — —  and  6  i 
X  rms 

distortion  at  the  aperture, 
tion  due  to  beam  guality  an 
input  parameter  1DFLM1.  DE 
eguaticn  3. 19  and  ther  uses 


s  the  BUS  value  of  the  phase 
Irel  us  the  intensity  degrada- 
d  is  egual  to  1/(N)*.  N  is  the 
I HA V  first  evaluates  a2  using 
this  value  as  a  starting  point 
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to  ccapute  a  acre  accurate  a2  by  iteration  using  subroutine 
EQIHE1.  a2  is  used  tc  commute  the  beam  quality  M1F  and  the 
team  quality  Irel  as  in  equation  3.18  .  a2  is  adjusted,  and 
the  process  repeated  until  the  Irel  found  by  equation  3.18 
is  equal  tc  that  deteimined  by  1/(N)2. 


|  TA  ELE  VII 

j  DETBAV  Program  Variables  and  Definitions 


1 

I 

i 

I 

1 

i 

i 

I 

J 

J 

t 


Variable 


6 

rms 


A 

1/  N  2 


For  t  ran 

III - 

?AB££ 
WAV  EEQ 

It  EiC 


Name 

Ilia 

7  A  EEQ 
3 AVEEQ 


Definition 

Constant 

Phase 

distortion 

parameter 

1/  (IDF1MI)  2 


6.  K1EE2  and  BH 

Ihese  tvc  routines  calculate  the  beam  quality  M1F 
array,  lie  tlF  is  £Bef.  26] 

Mb  ("p )  =  exp  ^-k 2  [a2  -  C^(p)J^  (3.20] 

where  CK(P)  is  the  autocorrelation  function  of  the  phase  and 
is  the  phase  variance.  C  ^ ( p)  is  assumed  to-be  Gaussian, 
letting  C  ^  <  p)  •=  a2exp[^- ( p/L)  ,  where  L  is  the  phase 
correlation  length,  results  in  [Ref.  27] 
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Mb(p)  =  exp 


2  TT  6 


rms 


[- 


1-exp  -  ( o/L) 


(3.  21) 


iritis 

— —  is  the  wavelength  HAS  phase  distortion  and  is  ar 
input  parameter  (WAVEEQ).  1  ,dlso  a  user  input  (SCAIrC), 
will  default  to  1/5  tie  diameter  of  the  aperture  if  ret 
otherwise  specified. 


TAE1E  VIII 

M1FEC  and  BM  Exogram  Variables  and  Definitions 


V  ariafl s 

~2tt6  "*2 

rms 


Via 


k  (m)  1 


Fortran  name 
3u5 


VAHEv 

F(I) 

N4 


“Earn 

Definition 

VA5EC 

Beam  guality 
variance 

BQM1F 

Beam  quality 
Ml  F 

N  4 

#  of  MTF 
integration 
incremen  ts 

- 

l  2tt/A  ) 
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c 


I2C1E3 


ISC1R3  calculates  the  zenith  iso^iana tic  angle  for 
use  in  the  adaptive  optics  portion  of  the  program.  Ihe  ancle 
is  given  by  £Hef.  28] 


0o  - 


.  314 


2.91 

6 .88 


sec  0 


h  4. 

atm 


t 


h 


V, 


(-• *2) 


For  the  rear  zenith  case,  sec0=l  .  c2  is  the  refractive 

n 

index  structure  constant  and  is  calculated  by  subroutine 
CN23  as  function  of  altitude.  See  Fried  £Ref.  293]  fcr  a 
discussion  cf  the  abcve  angle. 


TABIP  IX 

ISOTEB  Program  Variables  and  Definitions 


Variables 
0O  (rad) 


C 2  (a)  3 

n 


0  (deg) 


Fortran  Name 
3uE 


ISOANG 


C  N  2 


Ha  rn 
ISCAHO 


N  3 


OMEGA 


Definition 

Zenith 

isoplanatic 

angxe 

Refractive  index 

structure 

constant 

Integration 
intervals  for 
turbulence 

Zenith  angle 


1C.  SCI  NT  and  SI  NILS 

If  full  adaptive  optics  are  used,  amplitude  scintil¬ 
lation  effects  will  act  -to  degrade  the  perf cr nance  of  the 


73 


conventional  phase  correcting  adaptive  optics  system.  This 
subroutine  computes  tie  log  amplitude  variap.ee  of  the  scin¬ 
tillations  ror  zenith  traps mission.  The  approximation  used 
is 


hfhativ  s/6 

o2  =. 56k  /  h  C2(h)dh 

z  /  n 


(3. 23) 


lor  the  eff  zenith  case 


q2  =  o2  sec ( 0) 
z  z 


13.  24) 


where  0  is  the  zenitt  angle  £Ref.  30]. 
ance  less  due  tc  scir tillat ion  is 


Ine  relative  iiraci 


AMP  LOSS  =  exp( -crsec  (  9) 


<3  .25) 


Ihe  correction  fer  tie  off  zenith  case  is  applied  in  suircu 
tine  SINUS  in  the  ancle  loop  cf  the  main  program  when  full 
adaptive  optics  are  specified.  This  loss  factor  has  teen 
shown  tc  le  limited  tc  approximately  .55,  therefore,  a 
default  value  is  used  in  the  main  program  in  the  case  where 
eguaticr  3.25  yields  a  value  less  than  0.5  £Hef.  31]. 
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TABLE  X 

SCIMIS  and  SCIN1  Program  variables  and  Definitions 


Variaf ie 

X  (nf^ 

C2  ( K)  2/3 
n 


Fortran  Name 
3uE  H  am 

CK 

CN2E 


a 


2 

x 


SIGX2  SIGX2 


Definition 

2tt/ A 

Refractive  index 

structure 

constant 

Log  amplitude 
variation 


Amp  less  TAill  TAMP 


Loss  caused  ty 
scintillation 


11.  CN2K 

CN2B  computes  tie  index  of  refraction  structure 
constant  as  a  function  of  altitude.  is  determined  usin 
Eufnagel's  model  with  an  added  term  to  include  the  near 
surface  effects.  [Ref.  32]  [Ref.  33] 

Cn  =  (5.94xlO-SJ)h‘"exp(I^j  (3.26) 

+  (2.7xl0-'6)exp^^-)  . 

+  ll0""lexP  (lij§) 

Here,  tie  units  of  h  are  meters. 

1 2  .  E HC IRE 

Ibis  subroutine  computes  Yura’s  [Ref.  34])  lateral 
coherence  length, Po  ,  in  terms  of  Fried's  [Ref.  35]  coher¬ 
ence  cianeter,  ro- 
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in 


Po  =  r  o/2 . 1 


(3. 27) 


Substituting  Fried’s  definition  for  r 


P  o  - 


2 . 1 


O  Q1  ,  /  atm 

km  k2/h  c„(h)dh 

Jht 


_3 


Vc 


(3.26) 


TABU  XI 

EHOTRB  Prociaa  Variables  and  Definitions 


Varial i e 


k  (it) 


-  i 


-2 


C2  (E) 
n 


p0  (meters) 


Fortran  Name 
Sun  ain 

CK 

CN2 


HHC 

N 


riECO 
N  3 


£L 

2tt/X 

Rexractive  i ndex 
str  uc  ture 
constant 

Coherence  length 

#  of  integration 
intervals  for 
turbulence 


13.  Ill  and  MTFJI1 

311  AND  MTFJI'I  function  algorithmically  the  sane  as 
Eh  and  Mil  EC-  Together,  they  compute  the  jitter  3TF  array. 
MIFJI1  is  called  by  JIT  for  each  radial  increment.  The 
jitter  Mir  is  given  iy 


M,(p)  =  exp/k  p  (2ap} 

J  V  8 


(3.25) 
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T 


where  2q  is  a  user  specified  input.  As  with  teas  guality, 
P 

this  quantity  is  applied  ac  the  aperture. 


TAEII  XII 

JI1  and  aiFJIT  irogram  Variables  and  Definitions 
j  Variable  Fortran  Name 

~  Ini  ~ 

k  <mf  1  CK 

K_.{^)  F  JI1MTF 


Definition 
2  u/ A 

Jitter  MIF 


14.  IFIIHS 


Ihis  subroutine  calculates  that  part  of  thernal 
blooming  net  dependent  on  the  slant  patn  of  the  beam.  Ihe 
equation  used  by  GUISiSVG  for  the  phase  distortion  due  tc 
thernal  blooming  is 


A*(x,y>=  ^  fX‘  I(x’'y,dx'] 

L  a  Po  Y  J  J 


atm 


a (h) exp  (- 


(-sec©  / 

k 


a(h)+a(h)  dh. 


dh 


(3.20) 


‘t  V0cos£  +  wh 

Note  that  tie  first  term  in  this  expression  is  invariant 
with  respect  to  team  yath  while  the  second  term  contains 
such  path  dependent  varxables  as  wind  ,  wind  due  tc  slew, 
and  extinction  coefficients.  Ihe  first  term  is  evaluated  by 
this  subroutine  and  the  remainder  evaluated  later  in  the 
program  inside  the  sheet  angle  iteration  loop  by  subroutine 
AV. 
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•Standard  values  use 
2-  7  2  x  1 C  4  ,  1 . 0 1  x  1  Cs  J/a3,  an 

screen,  I  (i,j)  ,  is  construe 
intecral.  Note  that  program 
half  flare  cf  the  aperture 
t  is  the  aperture  dianeter. 
in  sutx'cutire  PHVAE  for  thi 


d  for  n  -1,  rot  an^  Y 
d  1.4,  respecti veiy .  k  phase 
ted  by  iteration  or  the  intensity 
performs  this  integration  ever  a 
defined  by  -b<y<b  and  0<x<t  where 
The  results  are  adjusted  later 
s  method. 


TA Elf  XIII 


PEIPHS 

Program 

Variables 

and  Definitions 

Variable 

Fortran 

III - 

Name 
“3  ain 

Definition 

k  ( m)~  1 

CK 

- 

2  7T  /  A 

n  - 1 

0 

const. 

Eef racti ve 
index  term 
(2.  72x  1  0 4) 

Y 

const. 

— 

Eatio  of  spa  cif 
heats  {1.4} 

Po 

c  onst. 

— 

Atmos,  pressure 
(1.01x10s)  N  /  m2 

4>  (i/Y) 

PH  (i,j) 

PH  (i,j) 

Thermal  bloomin 
phase  array 

15.  EESPHS  and  CCEFF5 


Cnee  the  phase  screen  has  been  constructed  by 
subroutine  EETPHS,  subroutine  EESPHS  removes  the  mean,  tiit 
and  focus  curvature.  Ike  result  is  the  phase  aberration  due 
to  thermal  blooming  alone.  Zernike  polynomials  are  utilized 
for  expressing  these  phase  distortions  with  subroutine 
COEFFS  providing  the  ieguired  ccexf icients.  These  polyno¬ 
mials  are 


t 
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20  (x  ,y)  =  a0  (mean)  (3.51) 

21  (x,y)  =  ai  x  (x  tilt)  (3.32) 

22  <x,y)  =  a2  y  (y  tilt)  (3.53) 

23  (x,v)  =  a3  lx2*y2)+alf  (focus)  (3.54) 

lie  expansion  coefficients  axe  given  below  and  are  computed 
relative  to  a  uniform  aperture  weighting  function  W(x,y). 
K(x,y)  =  1  inside  the  aperture  and  zero  elsewhere.  Ihe  ante 
grals  have  been  multiplied  by  2  to  compensate  for  the  hall 
plane  integration  of  subroutine  PHRPHS. 

f0=  2a0  (J)  (x#y)«  (x,y)  dxdy  (3.55) 

fi=  2ai  <j>  (x,y)W  (x,y)  xdxdy  (3.56) 

f2=  2a2  $  (x,  y)W  (x,y)  ydxdy  (3.57) 

f  =  2  (a  (x2  +y2)  +a  )  (x,y)  W  (x,y)  dxdy  (3.38) 

3  3  if 

Ihe  final  phase  correction  is  given  by 

<f>(x/y)=<Mx,y)-fiZ1(x/y)-f2Z2(x/y)-fjZ3(x/y) 

(3.53) 

-f  0Z  o (x,y) 


Por  a  discussion  of  least  sguares  fitting  of  Zernike  pclyrc 
mials,  see  £Ref.  36]. 
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16  .  I H  V  A  E 


Tie  phase  variance  due  tc  thermal  blooming  is 
computed  hy  this  suhicutine  using  the  residual  phase  screen 
provided  by  BESPES.  Normalized  with  respect  to  the  aperture 
field,  tie  variance  is  given  by 

2 

a z  _  U  (x,y)  j)2  (x,y)  +  U  (x,y)  tj>  (x,y)  (2.  QO) 

U (x, y) 


TABLE  XIV 

PHVAfi,  BESPHS,  ana  CCEFS  Program  Variable  Definitions 


Variable 
4><X/ y) 


a  -  a  , 

0  4 

V  *3 


Fortran  Name 
— H  am 


5uF 

PH  (i,j) 

S1G2QQ 

Al-  A4 

PHUEAN 
PHII1X 
PHI17Y 
PHE  CCUS 


PK  (i,j) 
SIGSQO 


Definition 

Pnase  array 

Pnase  variance 

Zernihe  coeffs. 

The  expansion 
coefficients 
relative  to  a 
uniform  weighting 
function 


17.  AV 

AV  evaluates  tbe  path  dependent  part  (second  term; 
of  tie  thermal  blooming  equation  5.30  .  (see  section  PBIPES) 
a(h)  and  cr  ( fc)  are  the  absorption  coefficient  and  scattering 
coefficient,  respectively.  V0  is  the  user  specified  wind  and  to 
is  the  slew  rate.  The  denominator  represents  the  total  tran¬ 
sverse  wind  component  across  the  'beam.  0  is  the  zenith  angle 
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and  £  is  the  angle  cf  attack  cf  V0  with  respect  to  the  leak. 

V  has  teen  assumed  tc  be  opposite  in  direction  and  parallel 

o 

to  the  target  to  transmitter  ICS  slew  motion.  Note  also 
that  V  Q  will  be  applied  as  a  constant  the  entire  length  of 
the  team  path.  This  is  discussed  in  chapter  one. 


TABIF  XV 

AV  Program  Variables  and 

Definitions 

V  arianc e 

Fortran 

Ill - 

Name 

Hain 

i 

a(h)  +  a(h) 

ALP 

- 

Total  extinction 
coefficient 

a(h)  (km)  1 

ALS 

- 

Absorption  coexf. 

sec  (6) 

SEC CFG 

SECCMG 

secant  or  zenith 
an  gle 

V  0  (a/sec) 

VO 

VO 

'wind 

PHI I Cl 

PHIICT 

Slew  rate 

b  ^  (meters) 
atm 

HATMC 

HAIMC 

Height  of  tt  e 
atmosphere 

h  (aeters) 
t 

HT 

H IRA  NS 

Height  of  the 
transmitter 

IE.  EICCM 

The  purpose  cf  BLOOM  is  tc  provide  reason  atie 
results  for  thernal  llccming  degradation  when  the  Irel 
values  are  telow  appicximately  0.3.  In  this  region,  the 
exponential  Strehl  relation  predicts  unacceptabiiy  severe 
results.  Therefore,  when  the  phase  variance  is  greater  than 
1.2,  EICCM  computes  at  Irel  value  based  on  GUTSMTF  results. 
GUTSMTF  is  a  full  wave  optics  propagation  cede  using  FFTs. 
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Curve  lit  p  cly  nomials  were  developed  using  GU1SKTF  results 
fcr  Lets  a  uniform  and  a  truncated  Gaussian  aperture  distri¬ 
bution.  Ihe  resulting  polynomials  are 

Irel  —  —  ~  (3.41) 

-.08705  +  2.9148a  +  .1723a2 


1 

Irel  =  - 

1.2877  -  2.6491a  +  4.09603a2 


{  J  • 


42) 


II  the  profile  under  consideration  is  a  truncated  Gaussiar, 
i.e.  the  aperture  diameter  greater  than  tne  waist  diameter 
cf  the  team,  then  the  routine  uses  the  truncated  Gaussian 
polynomial.  Otherwise,  a  co mtination  of  the  two  are  used, 

Ireltb  =  [l- (d/b)  2J  lrelu  +  (d/b)  2Ireltg  (3-43) 

where  d  is  the  aperture  diameter  and  b  is  the  waist  diameter 
cf  the  team.  [Bef.  27] 
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TA  ELI  XVI 

EIOGH  Program  Variables  and  Definitions 


Variable 

Fortran 

III - 

Marne 

~Harn 

Definition 

a2 

S 

SIGSC 

Phase  variance 

Il6ltb 

T 

TE1CCM 

Iotal  blooming 
produced  Irel 

lrelu 

THU 

Curve  fit  Irel 
for  uniform 
aperture  dist. 

lreltg 

TfiG 

Curve  fit  Irel 
for  Gaussian 
aperture  dist. 

I 

i 

i 


IS.  tlF  AIM 


C3FA1H  computes  the  atmospheric  MTF  and  if  speci¬ 
fied,  applies  a  tilt  due  to  turbulence  correction,  ihe  KTF 
is  giver  by 

Cp)  =  exp 


-(P/D)  ^ 


^ 1-d  (  p/D)  ^  (D/po) 


V, 


(3.44) 


where  L  is  the  aperttre  diameter,  p0  is  the  coherence  diam¬ 
eter  and  d  is  defined  as  1-ADAP.  AD  AP  is  the  fractional 
residual  tilt  due  to  turbulence  and  is  a  user  specified 
parameter.  If  ADAP=0,  M  will  represent  the  fully  turbulence 


tilt  cc 

r  rec 

ted 

MTF.  If 

a  d/  p  =  i ,  a 

will  be 

the  uncorrected 

MTF  and 

may 

be 

w  ri  1 1  e  r 

as 

r 

V 

(3. 

•45) 

MtCp)  = 

exp  1  - (p/po 

>  ] 

A£AP  nay  be  any  value  between  0  and  1. 
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TABLE  XVII 

£3  If  AIM 

Pro  cram 

Variables 

and  Definitions 

Variable 

For  tran 
Sub - 

Name 

~flarn 

Def inition 

F 

F  3 

Atmospheric  MTF 

po(meters) 

RHO 

RHO 

Coherence  di  a. 

£  (meters) 

D 

DIA 

Aperture  dia . 

2  C.  1  h  r  F  ZL 

•I£££EL  computes  the  approximate  Irel  for  a  tilt 
corrected  system.  This  value  is  used  for  ccmpari  son  uatn 
the  Iiel  produced  by  a  prefect  adaptive  optics  compensated 
system.  The  purpose  cf  this  comparison  is  to  determine  if 
noise  will  degrade  tie  AO  system  to  an  extent  as  tc  make  a 
full  AC  system  undesirable.  If  this  is  the  case,  the  program 
will  apply  tilt  correction  only  to  the  team. 

2  1.  I NDSHO 

Subroutine  FSLSHO  calculates  the  coherence  diameter 
that  wculd  result  in  the  Irel  value  acnieved  by  a  perfect  AO 
system.  I  his  .is  done  by  successive  calls  to  subroutine 
IEBBII  and  iterating  py  .  This  new  P0  then  becomes  a  factor 
in  the  AC  compensation.  Specifically  it  will  be  used  by 
subrcrtine  MTFATM  tc  produce  the  atmospheric  MTF  of  the  AC 
corrected  system. 
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TEE  Ell  and 

FFDBHC 

TABLE  XVIII 

Program  Variables  and  Definitions 

Variable 

Fortran  Name 

SuE  Ham 

Definition 

V?) 

F2 

- 

Atmospheric  K1F 

Ka(P) 

F  1 

IFBMTF 

Aperture  MTF 

Irel 

EX  D 

EXD 

Irel  of  perfect 

AO  system 

22-  1 13  ISO 

The  call  to  HUSO  is  made  when  adaptive  optics  have 
not  been  user  specified  and  iscplanatic  calculations  have 
net  teen  inhibited.  Also,if  the  signal  to  noise  ratic  at  the 
AO  sensor  is  such  that  the  AO  system  will  provide  tilt 
correction  cnly,  TLTIfC  is  called  ,  again,  provided  iscpla¬ 
natic  calculations  have  not  teen  inhibited. 

The  purpose  cf  T1TISC  is  to  include  beam  wander  due 
to  iscplana  tism.  The  2-sig  iia-p  tilt  is  computed  and 
combined  with  the  jitter  2-siyma-p.  This  combined  term  is 
then  used  tc  compute  the  jitter  MTF. 

22.  F El TOT 

FE1TCT  simply  averages  the  relative  intensities  due 
to  a  multiplicative  ard  an  RSS  (root  sum  squared)  approach  to 
combining  thermal  blccming  and  the  other  propagation 
effects. 


85 


t  n  Irel  +  Irel 
Irel  =  rss  m 


(3.^o) 


Irel  is  the  result  ci  tie  multiplicative  approach 
m 


Irel  =  Irel,,  *  Irel 
m  tb  o 


(3.  <I7) 


where  Irel^is  the  Irel  due  tc  thermal  blooming  and  lreIQ  is 
the  Irel  due  tc  team  cuality,  jitter,  turbulence  ,  iscplar.a- 
tism  and  adaptive  optics  effects.  I he  HSS  approach  is  given 

ty 


Irel 


rss 


(3.^8) 


TABLE  XIX 

B BITOT 

Pregram 

Variables 

and  Definitions 

Variable 

Fortran 

III - 

Name 

am 

Definition 

Iielt 

T 

TTCIAL 

Total  Irel 

Irel rss 

TR 

TESS 

irel  by  ESS 
method 

irelm 

TM 

TfiUII 

Irel  by 

multiplxcati.  ve 
method 

EXE 

EXD 

Irel  of  perfect 

AO  system 

lrel0 

— 

TTRUE 

Irel  due  to 

a^l  other  effects 

Irel  , 
tb 

TE1CCM 

Irel  due  to 
thermal  blooming 
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24.  ISOPLA 


This  routine  calculates  the  MTF  that  characterizes 
the  isoplanatic  effect  on  the  predictive  or  " look -ah ead " 
adaptive  optics  system.  Ihe  look-ahead  angle  is  the  major 
input  parameter  to  this  routine.  The  Fortran  code  was 
written  by  D.L.  Fried  [Eef.  38].  Fried  developes  the 
isoplanatic  dependency  of  the  AO  system  in  terms  of  the 
effective  antenna  gain  of  the  laser  transmitter.  The  MTF 
that  Fried  formulates  is  given  by 

t r 

r  -  z  -  5/ 

M.  (p)  =/  144.88-A-p  i  •sec(9)-y(vsec(0/p),<i>)dci>  (3.49) 

iso  / 


where 


x  =  vsec(0/p) 


(3.50) 


and 


Ic 


U  (x,  <{>)  =  /C^|l+(xh)  -%  1 1+2  (xh)  cos(f)+  (xh) 


-4- 


2(xh)cos  +  (xh)  1  6(ih 


■] 

<]i 


(3.51) 


0  is  the  zenith  angle  and  v  is  the  target  lead  angle. 

For  detailed  discussion  of  the  theory  and  explicit  develop¬ 
ment  of  the  Fortran  code,  see  [Ref.  39]. 
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JO  commutes  tie  zero  order  Bessel  function  based  cr 
the  input  argument.  :lhis  routine  is  called  oy  FAEFLC  in  the 
calculation  of  the  fai-field  ixradiance. 
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APPENDIX  A 


GUTSAVG  INPUT  FILE 


INPUT  DATA  HIE: 

LASER; 

CC  CSi  EDL  5P  ( 9)  THAHSITIOS 

4.9S21C  MICRONS 

CLIMATE: 

Blt-LATITODE  SU3HEE,  CLEAR  DAI 

CASE  : 

1 

DI  A 

IEIESCCEE  DIAMETEB 

HETEES 

0.  15C000E+ 01 

DIAOES 

CENTRAL  OB  SC  UBATICN  DIAMETEB 

METE  FS 

0.  15000  0  E+00 

BEAMSZ 

G  AC  SSI  AN  WAIST  TREU  TELESCCPE 

METERS 

C-  10  00  0  0  E+  02 

WAVE 

CAVIIY  WAVELENGTH 

METERS 

0.  4S921  0E-05 

PTOT AI 

AEESTCRE  TOTAL  PCWEfi 

WATTS 

C.  20  COO  0  E  +  07 

THSEE 

TCEECIEKCE  SEEING 

ARC-SEC 

0.0 

HGBNE 

EEIGE1  CF  GROUND  AECVE  MSI 

MET  c ES 

0 . 30  000  0  E+  03 

TDELKT 

TIMES  CIPPRACTION  LIMITED 

- 

0.  12 00 0 0  E+  0  1 

WA VEEC 

EES  A  V  ES  DISTORTION 

— 

C.  C 

SC ALEC 

PEASE  CORRELATION  LEMGTB 

METERS 

0.  30000  0  E  + 00 

HTRAhS 

I  BA  N  EM  1TTER  fcEIGET  ABOV E  MSL 

METERS 

C.3J0000E+03 

HS  AT 

S AIEIIIIE  AL1ITQCE 

METEBS 

C.  100000  E+  07 

THETKX 

M AXIMUf!  ZENITH  A  NGIE 

DEGREES 

C.  20  GO  0  0  E+  02 

LOFF 

EIJGE1  PATH  CFFSET 

METEBS 

0.  0 

SHOO 

YUBA'S  IURB.  COHERENCE  CIA 

METERS 

C.  C 

VO 

SIND  VELOCITY 

M/SEC 

C. 6QC000E+01 

SIC JIT 

2-SIGH A-P  JITTEB 

RADIANS 

0. 5 000 00 £-05 

ADAP 

llkC.  EESID.  TUBE.  TILT 

- 

C.  10C000E+01 

AOBICfl 

EICCMING  COR  RECTICN  (0-1) 

- 

0.  100000E+01 

AVGSET 

SPOT  AVERAGINF  FACTCc  (0-1) 

- 

C.  100000Z+  01 

NFLAGC 

INCLUDE  CLOUDS  MCEEI 

— 

0 

NFLAGA 

USE  EUII  ZONAL  A-C  SYST EH 

- 

0 

NOISC 

INKIEIT  ISOPLANATIC  CALCS 

— 

1 

ABSLCZ 

ZENITH  TRANS  FOB  A-C  SENSOR 

- 

0. T5C000Z+00 

XJT 

TARGET  RADIANT  INTENSITY 

W/STEB 

0.  25000  0E  +  03 

B5IDTE 

A-C  SYSTEM  BANDWIETH 

HERTZ 

C. 50GOOGZ+O3 

HA 

KUMEER  CF  A-C  SYSTEM  ACTUATORS 

0.  102400E+04 

N  1 

NUMEEB  CF  ANGLE  IMEBVAIS 

30 

N2 

AESCRP1ION  INTEGRATION  INTEBVA 

LS 

200 

H3 

TUBEULENCE  INTEGRATION  INTEHVALS 

100 

H4 

HIE  INTEGRATION  INTERVALS 

200 

H5 

THERMAL  BLOOWMING  I NTEGB ATICN 

INTER V  ALS 

100 
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APPENDIX  B 


GUTSAVG  OUTPUT  FILE 


LASER; 

CC  CB  EDI  5  P  (9 )  1EANSIT3CN 

4. 99210  MICRONS 

CLIMATE; 

MID-LATITUDE  SOMMER,  CIEA 

R  DAY 

CASE: 

1 

DIA 

TELESCOPE  E1AKETEP 

METERS 

.  150E  +  0  1 

DIAOES 

CSNIBA1  OBSCORAIICN  DIAMETER 

METERS 

• 15CE+00 

EEAMSZ 

GAOSSIAN  WAIST  THEO  TELESCCPZ 

METERS 

. 100E+02 

WAVE 

CAVITY  WAVELENGTH 

METERS 

. 499210E-05 

PTOT1L 

APZETUBE  TCTAL  FOWEH 

WA  ITS 

-  2  0  0  E  +  0  7 

THSEE 

IURBU1ENCE  SEEING 

ARC-SEC 

• 0CCE+00 

BGB  N  £ 

HEIGHT  OF  GROUND  ABOVE  MSL 

METERS 

•  3  00  E  +  0  3 

TDFLMT 

TIMES  DIFFRACTION  LIMITED 

- 

.  120E  +  0  1 

WAVEEC 

RMS  WAVES  £1 ST CRT  TON 

- 

- 00CE+00 

SCA 1  EC 

PHASE  COREEIATJCN  LENGTH 

METERS 

. 30CE+00 

HT  fi  A  fi  5 

TRANSMITTER  HEIGHT  ABOVE  MSI 

METERS 

.  3  00  E  ♦  0  3 

HS  A 1 

SATELLITE  A ITITOLE 

METERS 

•  1 OOE  +  C7 

THETBX 

MAXIMUM  ZENITH  ANGLE 

DEGREES 

.  3  00  E  +  0  2 

LCFE 

FLIGHT  PATH  CEISET 

MET  ERS 

.  00CE  +  00 

fiHOO 

YURA'S  TUBE.  COHERENCE  DIA 

METERS 

. OOCE+OO 

VO 

WIND  VELOCITY 

M  /SEC 

RADIA  NS 

. 060E+02 

SIG JIT 

2-SIGKi-P  JITTER 

. SC0E-C5 

A  LA  P 

FH  AC-  RESIL.  TORE-  TILT 

- 

.  1C0E  +  01 

AOB1CB 

BLOCM1NG  CCERICIICN  (0-1) 

- 

. 100E+C1 

AVGSFT 

SPOT  AVEBACINE  EACTOH  (0-1) 

— 

.  1 00  E  +0  1 

NFL AGA 

USE  FOIL  2CNAI  A-C  SYSTEM 

- 

0 

NOISC 

INHIBIT  ISC  PL AN  ATIC  CALCS 

— 

1 

ABSLCZ 

ZENITH  THANE  ECR  A-0  SENSCB 

- 

. 75CE+CC 

XJT 

TARGET  HALT  ANT  INTENSITY 

W/STER 

• 250E+03 

BSID1B 

A— 0  SYSTEM  EANLWILTH 

HERTZ 

.  5 0 CE  +  0 3 

NA 

NUMBER  OF  A-C  SYSTEM  ACTU  A  TCBS 

1024. 

N1 

NUMBER  OF  ANGIE  INTERVALS 

30 

N2 

ABSORPTION  INTEGRATION  INTIFVA LS 

200 

N3 

TUREUIENCE  INTEGRATION  INTEEVALS 

100 

N4 

MTF  INTEGRATION  INTERVALS 

200 

N5 

TREBMAI  BICCWMING  INTZGRATJCN 

INTERV ALS 

100 
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PATH  ANALYSIS  BESOLTS: 


STEF 

BANGE 

SLEW 

CMEGA 

T  IME 

AIM03.  TRANSMISSION 

fi  HO 

NC. 

(KM) 

(ME  AD/S  EC) 

(LEG) 

(SEC) 

SCAT 

ABSC5P 

(CM) 

1. 

999.738 

7.1582 

1.08 

2.  63 

0. 99802 

0.533  14 

51.  1 

2. 

1CC0.045 

7-1541 

2.16 

5.27 

0.99802 

0. 533C2 

51.0 

3. 

10C0.659 

7.1459 

3.24 

7.  SO 

0.99802 

0.53276 

51.0 

4. 

10  C  1. 57  9 

7.1337 

4.31 

10.53 

0.99802 

0.53242 

51.0 

5. 

1 0  C  2.  8  0  5 

7.1175 

5.39 

13.  16 

0.99801 

0.53195 

50.9 

6. 

10C4.334 

7.0973 

6.  4o 

15.80 

0.99801 

0.53135 

50.9 

7. 

1 0  C  6.  167 

7.0732 

7.52 

16.43 

0.99800 

0.53064 

50.  8 

8. 

10C6.301 

7.0453 

6-59 

21.06 

0. £9800 

0.52962 

50.8 

9. 

10  10.735 

7.0133 

9.65 

23.70 

0.99799 

0. 52888 

50.7 

10. 

1013.465 

6.97  8  6 

1C. 70 

26.33 

0.  597  99 

0. 52782 

50.6 

11. 

10  16.490 

6.94  0  0 

1  1.74 

28.  96 

0.99798 

0. 52665 

50.5 

12. 

1019-806 

6-8981 

12.79 

31. 59 

0. 99797 

0.52537 

50.4 

13. 

1G23-  4  13 

6.85  3  0 

13.82 

34.23 

0.  99796 

0. 52398 

50.2 

14. 

1027.305 

6.8048 

14.85 

36.  86 

0.  99796 

0.52248 

50.  1 

15. 

10  3  1-  4  e  1 

6.7537 

15.86 

39.49 

0.99795 

0.52088 

50.0 

16. 

1035.935 

6.6999 

16.87 

42.  13 

0.99794 

0.51517 

49.  8 

17. 

1040.  666 

6.64  3  5 

17.88 

44.76 

0.99792 

0.51735 

49.6 

18. 

1045.668 

6.5847 

18.87 

47.39 

0.99791 

0.51544 

49.5 

19. 

1050.939 

6.5237 

IS. 85 

50.  02 

0.99790 

0.51342 

49.3 

20. 

1056-474 

6.4606 

20.83 

52.  66 

0. 99789 

0.51  130 

49.1 

21. 

1062.269 

6.3955 

21.79 

55-29 

0.99787 

0. 50909 

48.9 

22. 

1068.319 

6.3288 

22.75 

57.  92 

0.  99786 

0.50679 

46.7 

23. 

1074.62  1 

6.2604 

23.69 

60.56 

0.99785 

0.50439 

48.5 

24. 

1081.169 

6.1907 

24.63 

63.  19 

0.99783 

0.50191 

49.3 

25. 

1067.959 

6.1196 

25.55 

65.  82 

0.99781 

0.49933 

48.  1 

26. 

10S4- 988 

6  .  04  7  5 

26.46 

68.45 

0.  99780 

0.49668 

47.9 

27. 

1102. 25C 

5.9744 

27.36 

71.  09 

0-  S  97  78 

0.49393 

47.7 

28. 

1109.740 

5.9006 

28.25 

73.72 

0.  99776 

0.49111 

47.4 

29. 

1  1  17.455 

5.8260 

2S.  13 

76.  35 

0. 99774 

0.48821 

47.2 

30. 

1  125.386 

5.7510 

5C.00 

78.99 

0. 99772 

0.48524 

47.0 
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THAKSEISSICM  ANA1YSIS  OUTIl’T: 


STEP 

LEANS 
DUE  TO 
SPREADING 

MISSION  COEEFICIEKTS 

AM  E  LOSS  TEERMAL 

WITH  AO  EICCilING 

RANGE 

SCALE 

MAX  IRRAD 
.  (KW/CM2) 

FLUENCE 

(KJ/CM2) 

1. 

0.  11826 

1 . 00  0  0  0 

0.01522 

0. 9999 

0. 5779E-04 

C.  1522E-03 

2. 

0.  1  1824 

1  .OOOOC 

C.O  152  1 

0.  9993 

0.57705-04 

C.  304  1  E- 03 

3- 

0.  1  182  1 

1.00000 

C.O  151  8 

0. 9981 

0.57535-C4 

C. 45555-03 

4. 

0.  1  1817 

1 . 00  0  0  c 

0. 01515 

0.9963 

0.57265-04 

C.  606  3E-03 

5^ 

C.  11810 

1.00  0  00 

C.0151  1 

0. 9938 

C .  5  692  E- 0  4 

C. 9562E-03 

6. 

0.  1  1802 

1.  OOOOC 

C.O  1505 

0. 9903 

0.5649E-04 

C.  S04  9E-03 

7. 

0.  1 1793 

1.00000 

C.  01495 

0. 9872 

0. 5  5  972-04 

C.  105 2 E-  02 

8. 

0.  1  1782 

1.  OOOOC 

C.O  149  1 

0.  9830 

0.5539E-04 

C. 1193E-02 

5. 

0.  1  1770 

1  .OOOOC 

C.O  1433 

0. 9783 

0.54725-C4 

C.  1 34  2 E- 02 

10. 

0.  1  1756 

1  .  00  0  0  c 

C. 01473 

0.9730 

0.5399E-04 

0.  1484E-02 

11. 

0.  1  1740 

1 . 00  0  0  0 

C.O  1463 

0.  9672 

0. 5320E-04 

C. 1624E-02 

12. 

0.  1  1724 

1  .OOOOC 

C.O  1452 

0.9610 

0. 5234E-04 

C.  176  2E-02 

13. 

0.  1  1705 

1 .00000 

C. 0  144  0 

0.9542 

0.5  142E-04 

C.  1898E-02 

14. 

0.  1  1686 

1.  OOOOC 

C-01427 

0.9470 

0.5  046  E-0  4 

C .  203  0  E-02 

15. 

0.  1  1664 

1 .00000 

C.O  141  4 

0. 9393 

0.  4  945  E- C  4 

C. 21611-02 

16. 

0.  1  1642 

1. OOOOC 

C. 0140  C 

0.93  13 

0.4839E-04 

C.  228  8  E-02 

17. 

0.11618 

1.00000 

C. 01385 

0. 9228 

0.  4  73  0  E-  0  4 

C. 2413E-02 

18. 

0.  1  1593 

1 .00000 

C.0136  S 

0.9140 

0.4618  E-0  4 

0.  253  4  E-02 

19. 

0.  1  1566 

1 .00  00  0 

C. 0 135 4 

0. 9049 

0. 4504E-04 

C.  26  53E-02 

20. 

0.  1  1538 

1 .00  0  00 

C. 01337 

0.8954 

0.4  387  E-0  4 

0.  27  6  8  E-0  2 

21. 

0.  1  1509 

1 . 00  0  0  0 

C. 01320 

0. 8857 

0.426  BE- 0  4 

C.  288  IE-02 

22. 

0.  1  1479 

1.00  0  00 

0.01303 

0.8757 

0.4 148E-04 

0.  29  9  0  E-02 

23. 

0.  11448 

1  .00000, 

C-0  1286 

0. 8654 

0 -  4  02  8E- 0  4 

C.  30  96  E-02 

24. 

0. 11415 

1 . 00  0  0  0 

0.0  1268 

0.  8550 

0 . 3  907  E-0 4 

C. 3 1 9  9  E-02 

25. 

0.  1  1381 

1-00000 

C. 01250 

0. 8443 

0. 3  786E- 04 

C.3298E-02 

26. 

0.  1  1346 

1 .00  0  00 

C. 01232 

0.8335 

0.3665E-04 

C.  33  95  E-02 

27. 

0.  11310 

1.00000 

C.O  121 3 

0.  8226 

0-3545E-04 

C.  34  8  8  E-02 

28. 

0.  1  1273 

1  .OOOOC 

C. 01195 

0.81 15 

0.3  427  E-0  4 

0.  257  8  E-02 

25. 

0.  11235 

1 .00000 

0.01176 

0. 8003 

0.3309E-04 

0. 3666  E-02 

30. 

0.11196 

1 . 00  0  0  c 

C. 01158 

0.7891 

0.3  193E-04 

C. 3750  E-02 
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CASE  LATA  AN E  CALCULATED 

FACTORS  : 

TRANSMITTER  DIAMETER 

- 

1 50. 00 

CM 

CAVIT’i  WAVJLENGTH 

= 

4.  9  92  1  OP 

OM 

APERTURE  TCTA1  POWER 

r 

2.0000 

MW 

TIMES  DIFERACTIGN  LIMIT 

= 

1.20 

{WIDE  ANGIE  SCATTERING) 
TRAN  SMUT  ER  ATTITUDE 

s 

3  00.  CO 

M 

SATELLITE  AITITUCE 

= 

10C0-  CO 

KM 

MAX  ZiNITE  ANGIE 

= 

30.  00 

DEG 

FLIGHT  PATE  CEfSET 

= 

0.  0 

KM 

RHOO 

= 

0.5105E+C2 

CM 

WIND  VELOCITY  • 

= 

6.  CO 

M/SEC 

HEIGHT  OF  CECDNC 

= 

3  00.  CO 

M 

OPT  SEEING  AT  55C0A 

= 

0.  1 5E  +  0 1 

ARC-SEC 

BL COMING  AIAPTIVE  OPTICS 

I ACTOR  =  0. 

1 OOOE+O 1 

IHEAL1ANC  E  AREA  AVERAGING 

fACTOR  =  0. 

1000E+0  1 

TELESCOPE  Cl M  ENSIGNS : 
OUTER  DIAMETER 

150.00 

CM 

INNER  DIAMETER 

= 

1  5 . 0  C 

CM 

GAUSSIAN  WAIST  DIAMETER 

= 

1000.00 

CM 

EH  S  SAVES  IISTCRTION  =  0.1031E+00 

PHASE  COEEEL1 ATICN  LENGTH  =  0.30  00E+02  CM 


FOB  TEE  F  C II  20NAL  A-C  HOCIL: 
ZENITH  TRANE.  AT  A-C  SENSOR 
TARGET  EAIIANT  INTENSIT 1 
ADAPTIVE  OPTICS  BAN  DWIDT  E 
NUMEER  Cl  ACTUATORS 


=  C . 7500  E+OQ 
=  0. 2  500  E+  03  W/STER 
=  0 . 5000  E  +  0 3  HERTZ 
=  1024. 


THIS  EUN  - 

IS  THE  EASIC  CCDE  WITH  CKI7  SOME  MEASURE  OF  TILT  CORRECTION 
AND  WITHCOT  AN  ISOPLANATIC  MODEL 


NOMBEE  OF  INTERVAL  STEPS  ECR: 
ANGULAR  INTERVAL 
ABSCRPT I C  N  INTEGRATION 
RHC  CA LCD! ATICN 
MTI  CALCULATION 
THERMAL  ELCCMING 

ZENITH  LOG  AMPLITUDE  VARIANCE  = 
RELATIVE  IEEACIANCE  REDUCTICN  = 

ZENITH  LOCK- AREAL  ANGLE 

ISOPLANATIC  JITTER  (2-SIGMA-P)  = 

2-sigma-p  Ream  jitter 

TURBULENCE  JITTER  REJECTICN 
(RESIDUAL  =  I NIT  AL  *  ADAP) 


=  30 

=  200 

=  100 

=  200 

=  100 

0. 4905E-02 
0.  995  1E+00 

0.5227E+02  URAD 

0. 2004E+07 

5.00  URAD 
0.  100QE+0 1 


========  PROPAGATION  RESULTS 

INTEGRATED  FLUX  CN  TARGET 
TOTAL  ILLUMINATION  TIME 


0.007499  KJ/CM2 
157.97  SEC 
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APPENDIX  C 


GUTSAVG  PROGRAM  LISTING 


c  *** *  ******* 
c 

GDTSAVG 


[***  *******  ************  : 


[************  i 


NAVAL  PCSTGEAEUATI  SCBOOI  VERSION  1.2 
MODIFICATION  1,2 
LA IE  SI  CHANGE  DATE  22  PEE  84 

THIS  VERSION  OF  GL1SAVG  HAS  EE E N  fCEIFIED  FOR  IBM  370/360 
CCMPAT ABI LITI.  VARIABLE  NAME  LENG1HS  HA VE  EZEN  CHANGED  AND 
A  F  PBOAPRI AT  I  FORTRAN  CHANGES  MADE-  MACHINE  GENERATED  ERRORS 
STILL  OCCUR  COE  TC  UNRESOLVED  CDC/I3M  DIFFERENCES 


C 

c 
c 
c 
c 
c 
c 
c 
c 
c 
c 

C  THE  CLOUD  MCEZL  BAS  EEEN  REMOVED  PROM  THIS  VZRSIO! 

c 

c 

Q ************************** ***************: 

c 

C== ========= =  =  =  MEMO  FI  ALLOCATION  AM  VARIAELE 

C 

REAL  BQKTF  (3  0  C  ) 

REAL  J IT  KTF  j  3  C  C) 

SEAL  A  (  ICC, 14) 

BEAL  B  (1 5) 

REAL  PHMC1.2C1) 

REAL  TI5C  <30  0) 

REAL  LOFE.MU 
REAL  IRRM1F  (30C) 

REAL  LY.IZ 

BEAL  NAbCET,NG,ISCANG,ISOANO, NA 
C 

INTEGER  TITLE1  {60) - TITLE2  (80)  ,TIT1E3  (80) 

INTEGER  EESCR  (4C,6Q) 

C 

COMMON  / A TMO/  EAT  MO 

COMMON  /EC/  TCFLM1, VARBQ, BAVEEC ,SCALEQ 
COMMON  /AR3FLE/  E I  A , LIAO ES,  EE AMSZ ,U0 


MAY 


<******** 


ASSIGNMENTS  ========= 


C==  === ======  =  ==  CONSTANTS 

C 

ML  =  20 1 
JJ=0 

T  5TOLD  =  0 • 

F  4=1 

PI=2,*ARSIN  (1 • ) 

MLE=ML/ 2 ♦ 1 
B  A  TMO=  3  .  E  4 
T  AU=0 . 

X  M  B= 1  - E+6 
XCM=  1.  E  +  2 
XKB= 1 . E-2 
XK2=  1.  E-7 
RTD  =  90  ./AESIN  (1.) 

R  E ARTH=  6 • 4E6 
MU=3.  9  86E 14 

1BETSF=2.*PI/ (24.*3600. ) 
F=BEABTH*TEET£F 


ASSIGNMENTS 


C__ 

c 


VSURF= 

C AIL  ERRSET  (208,256  ,-1,1,1) 

=========  DATA  INPUT  SECTOR  ======== =============================c 


READ 

READ 

READ 

READ 

READ 

READ 

READ 


(5, 510) 
5, 510) 
5,  510 
5,  5  10 
5,  510 
5,  5  10) 
5,510 


1TI TIE  1  (I)  ,1  =  1,80) 
TITIE2  (I  ,1=1,80} 
TITIE3  I)  ,1=1,80) 


CESCR 

EESCR 

CESCR 


CESCR 
i EESCR 
(CESCR 


.1, 

,11] 

1  ,H  = 

1,50) 

2 , 

,n 

1  ,H  = 

1,50 

3- 

,11 

1  ,H  = 

1,50 

r  II 

1  ,H  = 

1,50} 

5 

,n 

)  ,11  = 

1,50 

i 

rll 

1  ,H  = 

1,50 

(7, 

rH 

)  ,H  = 

1,50 

CIA 

DIACBS 
EEAKSZ 
BA  VE 
FTC  T  AL 
THSEE 
RGRND 


(EITHER  TEFLMT  CF  HAVEBQ  MUST  EE  DEFINEC,  BUT  NOT  BOTH) 
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nnn  oonnnnnnnon  oonnonnn  oo 


READ  | 

f  5 ,  510' 

LESCR 

;e,ii) 

,11=  1,50) , 

f 

TDFIMT 

READ  l 
READ 

5,5  10 
5,510 

LESCR 

LESCS 

#11=  1  /5o( 
f ,  II  =  1 , 5D 

I  t 

SAVE3Q 

SCALBQ 

READ 

5,5  10 

i  i 

IE  SCR 

11,11, 

,11=1,50 

1  9 

HTRANS 

READ 

5,  510 

LtSCE ! 

!  12,  II 

j ,11=1,50 

)  / 

HSAT 

BEAD 

5,510 

LESCR 

.13,11 

,11  =  1/50 

I  f 

THETMX 

READ 

5,  510 

LESCR 

14,  II, 

, 11=1,50 

|  9 

LCFF 

READ 

5, 510 

LESCR  ( 

15,11 

,11=1/50 

I  9 

RK  CO 

BEAD 

5,510 

LESCR  1 

.16,11 

,11=1,50 

)  § 

VO 

READ 

5,  510] 

LESCR  { 

17,  II 

i/II  =  1, 501 

F  / 

SIGJIT 

READ 

5,510 

LESCR 

'18, IT 

,11=1,50 

/  9 

ADAP 

READ 

5,510 

1  i 

LESCR 

19, II] 

,11=1, 50 

)  9 

AOBLOK 

READ 

5,  510 

LESCR 

20,11 

,11  =  1,50' 

1  9 

AVGSPT 

READ 

5,  520 

LISCB 

21,  II 

1/11=1,50 

)  9 

NFIAGA 

E  E  AD 

5,520 

LESCR 

22,11' 

l  ,n  =  i,5o; 

)  9 

NOISO 

IF  FULL  AC  CORRECTION  IS  i 

NOT  ESED, 

T 

HE  folic 
BWIDTH 

PARAMETERS  ARE  NCI  REALLY 

US  El -  ONLY 

HAS  AN  EFFECT,  it  CACSES  the  ICCK  AHEAD  ANGLE  TO 
BE  LARGER  BY  TAU*IRILOT. 


BEAD 

[5,510) 

(EESCE (23,  II) 

READ 

5,510 

(IESCE  24, II) 

READ 

5,510) 

(CESCB  (25,11) 

READ 

5  510 

CESCE  26, II) 

ABSLOZ 

XJT 

B  'WIDTH 
NA 


ITERATION  LOOP  UNITS: 

N 1  NUM3EF  OF  SDETEETA  INTERVALS  FROM  0  TO  THETMX 

N  2  NUMBER  OF  INTEGRATION  IN  T  E  E  V  ALS  FOR  ABSORPTION 

K 3  NUMBER  OF  INTEGRATION  INTERNALS  FOR  TURBULENCE 

N  4  NUMBER  OF  INTEGRATION  INTERVALS  FCR  MTF 

N5  NUMBER  OF  SUEINTERVALS  USED  FOR  SLANT  PATH  UPDATE  FCR 

THERMAL  ELOO  MING 


BEAD 

READ 

READ 

READ 

READ 


(5,  5  20 
5,520 
5,520 
5,  520 
5,520 


27,  II)  ,  11  =  1, 50)  ,  N1 

28.11  ,11=1,50  ,  N2 

29,  II)  ,11  =  1,50)  ,  N3 

30,  II)  ,11=1,50  ,  N4 

31. 11  ,11=1 , 50 i  ,  N5 


ECHO  CHECK  OUTPUT  SECTOR  ================== ======C 


WBITE 
W  BITE 
WBITE 
WRITE 
WRITE 
WRITE 
WRITE 
WRITE 
WRITE 
WHITE 
WRITE 
WRITE 
WRITE 
WRITE 
WRITE 
WRITE 
WRITE 
WRITE 
WRITE 
WRITE 
WRITE 
WRITE 
WRITE 
WRITE 
WRITE 
WRITE 
WRITE 
WRITE 


i#  ec) 
i  ,ec) 

I  ,sc 

11=1,50) , 
11=1,50  , 

II  =  1,50)  , 
11=1/50  , 
11=1,50  , 

_ _ 11=1,50  , 

7,  H  i  ,11  =  1,50  , 
~  11=1,50  , 


) ,11=1 

JU  1  , 

,  50 

,11=1 

,50 

/ 1 1  =  1 

,  50 

,11=1 

,  50  i 

,11=1 

,50 

,11=1 

,50) 

,11=1 

,  50  i 

,11=1 

,50 

,11=1 

,  50 

,11=1 

,  50  , 

,  11=1 

,  50  i 

,11=1 

,  501 

!  ,H=1 

,50 

,H=1 

,  50 

,H=1 

,50 

LI  A 

DIAOBS 
EEAMSZ 
WAVE 
PTOTAL 
IH  SEE 
HGBND 
TDFLMT 
WAVESQ 
SCALBQ 
HIRANS 
HSAT 
1 EET  MX 
,  LOFF 
,  RHCO 
VC 

S IG JIT 
,  ADAP 
,  ACBLOM 
AVGSPT 
,  NFIAGA 
NCISO 
,  ABSLOZ 
,  XJT 
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nnnnnn  nonnn  nnnnnn  nnnnn  nnnnn  nnnnn  nnnnn  r»  n9n  n<i?n 


3WIDTH 
N  A 
N  1 
N  2 
N3 
N  4 
N  5 

==  =  ==========  DATA  REALIGNMENT  SEC1CE  ======================== =======C 


WRITE 

[6 ,180] 

)  {DESCR  i 

;25  ,n; 

1 ,11=1,50) 

WRITE 

6,  130 

1  {CESCE 

26  ,  ii; 

1 ,11=1,50) 

WRITE 

‘6,185 

)  IDE  SC  H  J 

2*7 ,11 

I ,11=1 ,50) 

WRITE  { 

6 , 185] 

1  DESCR  ( 

28,  II 

I  ,11=1,50) 

WRITE 

6 ,185 

)  (DESCB  | 

29  ,11' 

I ,11=1 ,50) 

WRITE 

6,  185 

1  (DESCR  | 

3°,  II 

1 ,11=1,50 

WRITE 

6,185 

DESCR  | 

31  ,11 

,11=1 ,50 

(SCA1EQ.SQ. C. )  SC  AL3 Q= El  A/  5  . 

(WA  VEEQ . N  E . 0.)  VAEBQ=  (2. *PI*WAVEB C) * 
J±TO=SIGJIT 


If 
If 

SGJ1  _ 

If  (NFLAGA.  2£  .  1)  ACAF=0. 

IF  f  NOI S  C . EQ .  1 )  ISCAK0=1-E10 
MASQRT=SCET]N A) 

LOGA3S=-3lOG  (AESLC2) 

110=1. /EWIDTH/PI 

==============  INITIAL  CALCULATIONS 

IF  (RHGO.EQ. 0. AND.IHSZE.  HE. 0. ) 

1^HETMX  =  T?:fTMX/5TD 
HSA1D=HSAT-HTB ANS 
HAEVGD=HT5ANS-fcGEND 
IF  {HABVGL.LT,  C.  )  H  AEVG  D=0. 

R=REARTB*ETRANS 
RS=REABTH*HSAT 
NFIGSN  =0 


*2 


BHOO  =  (WA VE/. 55E-6)  **  1.2* ( -  0  54/THSEE 


COMPUTE  EARTH  CENTER  ANGLE  OF  E-SET 

ANGOFF=LCFF/R 
A  2= ANGO  E  E/2- 

COMPUTE  COORDINATE  TEANSLATIONS  DUE  TO  CFF-SET 


LY=+2.*R*SIN{Al)*COS(A2) 
12  =  2-*  R*SIN  {  A  2 )  **2 


COMPUTE  EARTH  CENTER  ANGULAR  RATE 


TBETDT=SCBT (MU/BS**3) 


COMPUTE  CEEITAI  SFEEE 


V  SAT=THETDT*RS 


ADJUST  INTENT A  NEC US  SLEHRATE  ECR  ROTATION  OF  EARTH, 
assuming  both  are  cclinear. 


THETDT= AES (T H E 1 LT-T fi ETS F) 


COMPUTE  VELOCITY  CF  SATELLITE  RELATIVE  TC  THE  TRANSMITTER  SITE 


VS=THETDT*RS 


COMPUTE  INITIAI  LIMITS  ON  THE  TIME  AND  ANGLE  EARTH  CENTER 
.for  an  or  line  flight  path 
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EC  At!  G=  A  R  CCS  ( fi  *  S  I  N  (I  H  ETtt  X )  /R  S)  +  THETtiX  -PI/ 2  . 
TIMEX2-=EC  ANG/T  EZ1ET 
TIHTOT=2. *TIMHX2 
DTIM£= TI £EX2/ N  1 


IF  SATELITE  A1I1 ITUDE  IS  AS  TEE  ORDEB  Of  THE  RADIUS  OF 
THE  earth ,  gut sip  (footprint)  WORKS  SETTER  BY  SETTING 
the  time  step  instead  ox  the  aigle  step 

IF  (HSAT.GT.  REABTE/2.)  DTIME=  10. 

COMPUTE  VACUUM  IBEAEIANCE  AMD  APERTURE  MTF .  NOTE  THAT  THE 
PCINT  SPACING  HERE  MOST  BE  THE  SAME  AS  USED  IN  THE  MTF 
INTEGRATION  3EICW.  HENCE  WE  EEFINE  DX  NOW. 

DX=DIA/N4 

THE  CALI  TO  UOCST  DETERMINES  TEE  CONSTANT  THAT  MAKES  THE  EXIT 
APERTURE  POWER  FTCTAL. 


CAII  UOCST  (PTCTAI,M1) 

CALL  FARFLD  (B X , T JSC ,IB R MTF , M I , N4 , HS ATD , W AVE , PTOT AL , PMAXO ) 


IF  THE  WAVE3Q  IS  NCT  SPECIFIEt  AS  INPUT,  COMPUTE  ON  BASIS  OF 
TDFLMT. 

IF  (WAVEEQ.EQ. C.)  CALL  DETWAV  (IERMT F  ,  DX  ,  N4) 

IF  TDFLMT  WAS  NCT  SPECIFIED,  COMPUTE  IT  . 

IF  (TDFLMT. NE.  C.)  GO  TO  10 
CALI  BQIEEL  ( I B E MI F , N4 , D X ,T BQ ) 

TEELST=SCET ( 1 -/TEC) 

CONTINUE 


DEFINE  BF AM  QUALITY  MTF  ARRAY. 
CAII  BM  (EQMTF  ,DX  ,N4 ) 


DEFINE  JITTER  MTF  ARRAY. 

C AIL  JIT  (JITHTF,EX,N4, SG  JIT  C  ,  W AVE) 


COMPUTE  ZENITE  AESORPTION,  SCATTERING,  AND  RHOO  VALUES. 


CALL  A  BS  CEB  (  K 2  ,  HTB A  NS, T ABSO ) 
CALL  SCAT  (N^ , ETHANS, TSCATO) 

IF  (BHOO.FQ.O.)  CALL  SHOT  SB  (NS, 


HTRA  NS,  EAEVGD,  WAVE,  RHOO) 


DEFINE  SEEING  AT  fSOOA 


THSEE=  (WAV E/ . 5  5  E  -  6 ) ♦♦I. 2*  (. 05^/ RHOO) 

IF  (NFLAGA.GE.  1.  AND- NOISO.EQ.  C)  CALL  ISCTRB  (N 3, HTR ANS , HAB VGD 
1  , ISOAN0) 


,  WAVE 


CALCULATE  LOG  AMPIITUDE  SCI  NT  II IATIO K  FCR  ZENITH  ANGLE. 
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CALL  SCIM  (SIGX2  ,N3  ,HTR  A  NS  ,HGEND,WAVE) 


COMPUTE  HE  THERMAL  ELOOM  IN  G  PEASE  DISTORTION  THAT  does  not  depend 
SLANT  PAIB  CHARACTERISTICS 


CALL  PRTFHS  (  HI  ,MLH , WAVE , FK) 

CALL  RESFES  (Ml  ,ML H  ,  PH) 

CALL  PHVAB  { M 1 , M LH , P TOT AL ,P H , SI GSQO ) 


INITIALIZE  ACC  UMUI AT  ED  FLUENCE  TC  2EEC 


E1UX=0 . 


SET  SNR  E LAG  TC  CEE 


NELGSN  =0 


************  MAIN  PROGEAM  ICOP  ***************** C 


L  COP  ON  IREADIANCES  AND  ACC  UM  Cl AT  ED  ELDENCES  FOE  A  FIXED  LOFE 


DO  120  1=1, N1 


COMPOTE  TCTAL  LENGTH  OF  ILLUMINAIION  TIME  FROM  (+,-) 

OVER  P  ASS  ANGLE 

THIS  IS  THE  TCTAL  T I  ME.  RE  OSE  HALF  OF  THIS  TO  DETERMINE  THE 

FUNCTION  EVALUATION  ANGLE. 


TIME=2 . *  (  (I-  1)  *  IT Ifl  E  +DT I M E/2 .  ) 


COMPUTE  EARTH  CENTER  ANGLE  AT  CN  LINE  COORDINATES  AND  TIM  E/2. 
THIS  ANGIE  IS  TC  TEE  MID-POINT  OF  THE  INTEGRATION  INTERVAL. 
THE  ANGLE  ECANGE, COMPUTED  L  AT  EE ,  IS  TO  THE  UPPER  LIMIT  OF  THE 
integration  liHit. 

******  THE  COMMENTED  ECANG  SHOULD  EE  USED  WEEN  IT  IS  DESIRED  THAT 
THE  INTEGRATION  SHOULD  START  AT  THE  LCfi  POINT  AND  INCREASE 
TOWARDS  THE  ZENITH.  THE  CURRENT  CALCULATION  STARTS  AT 
THE  ZENITH  ABC  GCES  DOWN. 


**  ECANG=  (TIBTOT-TIBE) /2.*THETDT 
ECANG=  TI ME/2  .  *  TE  ETDT 


COMPUTE  TARGET  COCEDINATES 


X0=BS*SIN  (ECANG) 

Z0=RS*COS  ECANG 
X  =  X0 

Y=IY*COS  (ANGOE  E)  *  (ZO-R+LZ)  *  S I  N  (ANGOFE) 
Z=-LY*SIN  (ANGCFF)  +  (ZO-R  +  LZ)  *CCS  (ANGOEF) 


COMPUTE  CN  LINE  CE  SIGHT  ANGLE  OF  SAT. 

THIS  ANGIE  IS  NCT  USED  AT  PRESENT  .  IT  IS  THE  ANGLE  CF  THE 
SAT  AS  MEASURED  E  EC  M  A  POINT  UNDER  THE  GROUND  TRACK  . 


TEETA=ATAF{RS*SIN  (ECANG)  /  (RS*CCS  (ECANG)  -R)  ) 
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c 

C  CCMPOTE  BAN GE  1C  1 A  E  GET 

C 

RANGE=SQ£I (X* ♦ 2* Y** 2*1* *  2) 

C 

C  INSTANTANEOUS  SLEW  BATE 

c - 

c 

VX=VS*COS  (ECANG) 

VY=-VS*SIN  (EGA  NG)  *SIN  (ANGOFF) 

VZ=-VS*SIN (ECANG)  ♦CCS (AN  GOFF) 

iX=Y*VZ-2*VY 

SY=Z*VX-X*VZ 

WZ=X*V Y- Y*VX 

c 

PHIDOT  =  SCBT (WX**2*W Y **2 + W Z* *2 )/BANGE * *2 
C 

C  COMPUTE  ANGLE  EfCE  ZENITH 

c 

C  OMEGA  EQUALS  ANGIE  DCWN  FROM  2ENITH 

C  OMWIND  EQUALS  ANGLE  CF  WIND  A11ACK  TC  LCS  I?  TARGET 

C  MOTION  AN!  WINE  ABE  COPLANAR. 

C 

OMEGA= ABCCSfZ/EANGE) 

CMHIND= ARSIN (I/BANGE) 


COSOMG=CCS(OMEGA) 

SECOMG=1  ./COSCMG 
CCSWND=CCS (OMWIND) 

C 

C  ADJUST  FCE  SLANT  BATH.  BLOOMING  LOSS  ,  AESCBPTION  LOSS,  SCATTERING 

C  LCSS,  AN E  TURECIENCE  RHO 

C 

c 

C 

CALL  AV  (N5/YC#BHIDC1,HTBANS, CCSWND, cosoog , E) 
c  SIGSQ=SIGSQO*E 

C  APPLY  THEBMAI  E1CCKING 

C  ADAPTIVE  OPTICS  DEGREE  OF  COMPENSATION 

C 

SIGSQ=SIGSQ*ACEICM 

C 

SIG=SQRT  (SIGSC) 

TELOOM=EXF f-SIGSQ) 

IF  (SIGSC.GT.  1.2)  CALL  BLOOM  (Z I A ,  BE  AMS  Z  ,  SIGSQ  ,T  BLOOM) 

TAES=TABSO**SECCMG 

TSCAT=TSCATO**SECCMG 

BHC=HHOO  *COSCMG**  (.6) 

C 

C  COMPUTE  LOOK  AHEAD  ANGLE  ASSUMING  EABTH  ROTATION  AND  SATELLITE  MOTION 
C  ARE  IN  THE  SAME  DIRECTION. 


TWCT=RANGE*2./3. EE 
DBT=TW  OT  *  (ABS  (BSAT-^SURF)  ) 
PRC JAN  =  A  ECOS  it/ EANGE) 

PEC J=DR1  *  SI N  (PBCJ AN) 

NO=PRO J/ RANGE +TAU*PH IDO T 
IF  (I.EQ.1)  ZENNU  =  NU 


C  IF  executed,  the  following  is  a  full  ao  simulation 

C 

IF  (NFLAGA.EQ . C)  GC  TC  40 
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c 

C  MODIFY  EEC  FCE  EfFECT  OF  AD AP 1 IV E  OPTICS 

C 

RHCLD=RHC 
B  BCU  =  R  BO 
R0=2.  1  *  R  BC 
C 

C  DETERMINE  APPROXIMATE  I-EEL  TC  SEE  IF  AC  SHOULD  BE  USED, 

C  given  PE ESENt_NGI£E . 

C 

DCE=DI A/EHO 

CALL  TRBEEL  (I  E  R  M TF  , LIA  , N  4,  D  X , FXDTST ,SHC) 

C 

C  COMPUTE  RESIDUAL  VARIANCE  DUE  10  PERFECT  ADAPTIVE  OPTICS, 

C  INFINITE  EANDfc'lITH,  AND  FINITE  NUMBER  OF  ACTUATORS. 

D  1  =  -32  0*  JEIA/H 0/N  ASQET)  **  (1  -  6  667) 

TES=LOGAES*SZCCflG 
TSENSR  =  E  Xr  (-1ES) 
c 

C  DEIERMINZ_PHA£H  V ARI£NCE_ ASSOCIATED  SITU  SENSOR _ _ _ _ 

PHSERR=8.5E-o*EWIDTH*NA**2* (RANG  E/5.  E5*2. 5/DI A) **4*  {  80 . /XJ T *. 5/ T SE 
1 N  SB) **  2 

C 

C  INCREASE  RESIDUAL  PEASE  VARIANCE  AND  COMPUTE  THE  I-REL. 

D  1  =  D1  +  PHSERE 
C 

C  COMPUTE  I -ESI  CF  AC  TURBULENCE  CORRECTED  BEAM. 

C 

EXD=EXP  (-D1) 

C 

C  IF  AO  CORRECTED  I-EEL  IS  LESS  THAN  NC  AC  CORRECTION ,  ASSUME 

C  ONLY  TILT  CORRECTION  IS  USED. 

C 

IF  (EXD-LT.EXD1SI.OB.NFLGSN.EC-1)  GO  TO  20 
C 

C  AIL  FBDEBO  {  E 2 D  ,  E HO , DI A , EHO U  ,  C X , N4 , IRE MTF) 

C 

GC  TO  30 
20  CONTINUE 

C 

C  THIS  SETS  A  FI  AG  TC  PRINT  THAI  AO  SYSTEM  IS  ONLY  TILT 

C 

RHC=RH  OLE 
NFLGSN— 1 
C 

30  CONTINUE 

C 

C  ESTABLISH  TEST  INCREMENTS  TO  SEE  IF  FULL  ISOPLANATIC  CALCULATION 

C  NEEDS  TO  BE  DC  N  E . 

C  COMPUTE  THE  ISCPLANATIC  ANGLE. 

C 

ISCANG=CCSOMG**  (2-6667)  *ISOAN0 
TSTNEW=NC/ISOANG 

TSTDIF  =A  E  S  (TSIKEN**!  .  66  66  67 -T  ST  CL  D) 
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c 

c  — 
C 

c — 
c 


40 

C 

c- 

C 

c- 

C 

c 

c~ 

c 

c 

c 

c 

c 

c- 

c 


50 


C 

60 

C 

c — 

C 

C 

c 

c  — 

c 


c 

c  — 
C 

c  — 
C 


c 


c  — 

c 

c 

c  — 


c 


70 


90 

100 


ETSDIF=EXF (IS  T  E  IF  ) 

IF  (NISO.EQ.  1  )  TSTOLE=TSTNEW**1.  666667 


DETERMINE  IF  ISCFLANA1IC  ST  F  SECULD  EE  RE-COMPUTED 


NISO=0 

IF  (ETSDIF.GT.  1.  1)  N1S0=1 
IF  (I.EQ.1)  N I SO=l 
if  ( noise. eq . 1 )  niso=0 
CONTINUE 


REDUCE  PEAK  IERAD1ANCE  BY  SLAM  DISTANCE 


T  MAXSC  = (ESATE/BAKGE)  **2 


IBIS  ALLOWS  ONE  VC  INCLUDE  TILT  I SO PL A  NAT  IS  M.  NOTE  HOWEVER ,  THAI. 
IF  THE  FULL  ADAFI1VI  CFTICS  ALGORITHM  IS  CSED  WITH  THE  ISO? L A N AILS K 
CORRECTION.  THEN  IBIS  CALL  SHCUIE  NOT  EE  US  ED- RE PE AT- S HOULD  NOT  BE 
USED.  ALSC  NOTE  TEAT  A  VEETICAL  TURBULENCE  PROFILE  IS  NEEDED 
WHEN  EVES  IBIS  FCUT1NE  XS  USED. 


IE  (NFLGSN.EQ.  1.AND.NOISO.EQ.C)  GOTO  50 
IF  (NOISC.EQ.  1-CE.NFIAGA.  EQ.  1)  GOTO  60 
CCNTIN  UE 
OMX=NU 
C  M  Y=0 . 

CALL  T  LT ISO  ( C I A , CSX , OS Y , BEST  II , SECO MG , HTR AN S , HAB VG D) 

SIGJIT=SQBT (SG JITC**2+RZSILT) 

CONTINUE 


COMPUTE  IRBADIANCE  DECREASE  DDE  TO  AMPLITUDE  SCI  NTH  LA  TI  ON  FOR  THE 
OFF  ZENITH  PASS.  NCTE  THAT  TEIS  AS  LOSS  HILL  ONLY  EE  APPLIED  IF 
FULL  ADAFTIVE  CFTICS  COMPENSATION  IS  ASSUMED. 


CALL  SINUS  (TAMF,SIGXZ,SECOHG) 


CALCULATE  EFFECTS  CF  JITTER  ANT  TURBULENCE 


SUM=0. 

DX  2=DX/2 . 


IT  IS  I MFCRTA  N  NT  TC  NOTE  THAT  ALL  THE  AfiBAYS  ARE  DEFINED 
TC  BEGIN  AT  DX/2. 


X=DX2 

DO  110  J=1,H4 
IF  (X.GT.DIA)  GO  TC  100 
F  1=IRR MTF  <J) 

CALL  MTF  ATM  (X  ,£I A , BHO, AD AP , F  2) 

F3=JIT  MT  F  (J) 

I?  (NFLAGA.EQ.  1 .  A ND .  NFL G S  N.  E Q  .  0}  GO  TO  70 
IF  { NO  ISC. EQ.  0)  CALL  MTF J IT  {  X ,SIG JI T, W A VE , F3) 

CONTINUE 
F  4=  1 . 

IF  J'NFLAGA.EQ.G.CE.  NCISO.  EQ.  1  .CR.  NFL  GS  N  .  EQ  -  1 )  GO  TO  90 
F4=TISO  ( J) 

IF  (NISO.EQ.  1)  CALL  ISOPLA  ( J  J  ,  NU  ,  X, SECOMG,WAVE, F4  , HTRANS, HAEVGC  ,T 

CCNf INUE^ 

F5  =  BQMTF  (J) 

SUM=SUM+E  1*F2*F3*F4*F5*X 
CONTINUE 
X  =  X  +DX 
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CONTINUE 

TIURB=SU  M*2.  *PI*CX 


1  10 


CALCULATE  THE  EES Cl I I NG  I NSTENTANEOU  5  INTENSITY 


CONVERT  TC  AN  APP50CE  WHICH  AVERAGES  THE  RESULT  OF  AN  RSS 
TREATMENT  OF  THERMAL  BLOOMING  WITH  A  MULTIPLICATIVE  TREATMENT 


SSQT=1 ./T1URE- 1. 

SSQB=1 - /TE10CM- 1. 

SSC-SSQT+SSQE 

TKULT=TTUS3*TS1CC* 

TE  SS- 1 ./  ( 1. +  S  SC } 

CALL  RELICT  ( T TOT AL , IRS S , TMU LI ) 
PMAX=PMAX0*TAES*1SC AT*TMAXSC*TT 


CT  AL 


APPLY  AVERAGING  EACTOR 
APPROPRIATE  TC  SPCT  SICE  LESIRED 


PMAX=PHAX*AVGSPI 


IF  FULL  AO  IS  DSED,  LECEEASE  IRSADIA NC£  FOB  AMPLITUDE 
SCINTILLATION  EFFECTS 


TAMPL= 1. 0 

IF  (NFLAGA.EQ.  1  .  A  KD  -  NFLG  S  N. 
IF  (TAMPL.LT-  0.5)  T  A  CPL  =  0.5 
P  MAX=P  MAX*TAMPL 


EQ.O) 


TAMPL=T AMP 


ACCUMULATED  FULX 


FLUX=FLUX+PMAX*CTIME 


COMPUTE  INTEGRATED  TIMES  AND  ACCUMULATED  ANGLES.  THE  IBP.ADIANCE 
FUNCTION  HAS  EEEN  EVALUATED  AT  INTERVAL  MID-POINTS. 


ECANGF= (TIME+E1IME) /2.*THETDT 
X0F=RS*S1N (ECA  KGr) 

Z0F=RS*CCS (ECANGF) 

XF=X0F 

YF=LY*COS  jANGCFF)  +  (ZOF-R  +  LZ)  *SIN  (ANGCFF) 
ZF=-LY  *S  I  A  (A NGCFF)+(ZOF-B+LZ)*CCS {AN  GOFF) 
BANGEF^SCET (XF**2+YF**2+ZF**2) 

OMEGAF=  AECOS  (ZF/RANGEF) 


TIME  OF  IILU  MINATICN  FROM  ZENITH  TO  TEETA-  THIS  IS  ALSO 
EVALUATED  AT  TEE  LOWER  LIMIT  AND  NOT  THE  MID-POINT. 

IF  THE  COMPUTATIONS  ARE  DONE  FROM  THF  THETMX  TO  ZENITH  INSTEAD 
OF  THE  WAY  THEY  ARE,  THE  TIME  DEFINITION  NEEDS  TO  3E  CHANGED. 


TlLLUM^TIfiE/2 . ♦CT1ME/2. 

=============  LOOP  RESULTS  STGRAGF  SECTOR  ========================== =c 


■1:21=5 

1/5  _ 

1.6  =T  SCAT 

1.7  =T  AES 


A 
A 
A 
A 
A 
A 
A 

A  jl',8)  =R  EC*XC K 


RANGE*  XKS 
=  P  FIDO!/ XK W 
=0  KEGA  F  *  BT  D 
=TILLUM 
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FINAL  RESULTS  OUTPUT  SECTOR 


E  (1)  =DIA*XCM 
E  (2) =WAVE*XM W 
E  3  =P1CTAL/XMW 
B  (4)  =TDF1  £T 
E  (5  =HTRANS 
E  6)  =BSAT*XKW 
E  (7  =THET«X*fiTE 
B  8  =L  OFF  *  XKW 
E  9]=BHOO*XCM 
B  it))  =  VO 
B  il  1)  =  HG  E  ND 
B 1 12)  =  IB  SEE 


WRITE  (6,220) 

WRITE  6,240)  (  (A  (I  ,  J)  , J=  1 , 8 )  ,  1=  1 ,  N  1  ) 

UCTTr 


WRITE 

WBITE 

WRITE 

WEITE 


El A  =  DI A*  X  CM 

DIA03S=DIA0BS*XCfi 

EEAMSZ=£EAM3Z*XCfl 

WRITE  (6,140)  EIA,EIAOBS, EE A MS 2 

SC ALBQ=SC A1BQ  MCE 
WRITE  (6,150)  WAVEEC  ,SCALEQ 
B  ( 1)  =  A  B  S I CZ 
B  (2)=X JT 
B  (3)  =B WIDTH 
B  (4)  =NA 

WHITE  (6,380)  (B  (II)  ,LL=  1  ,4) 

WRITE  6,260) 

IE  (NFLAGA-EQ.  0.  AKD.  NOISO.EQ.  1)  WEITE  (6,310) 
IF  (NFLAGA.EQ.  UAFC.  NCISO.EQ.O)  WRITE  (6,280) 
IF  (NF  LAGA-EQ-  1-ANC- NOI SO  .EQ .  1)  WRITE  6,290) 
IF  NFLAGA.EQ. 0.AND. KOI SO .EC . 0)  WRITE  6,300 
WHITE  (6,340)  * 1 , 1 2 , N3 , N 4  ,N5 

ESIGXZ-EXP (-SIGXZ) 

WEITE  (6,160)  SIGXZ , ZSIG  X  Z 

ZENNU=ZENNU*XM W 
WHITE  (6,320)  ZERKU 

BESTLT=SCBT (HESTLT) MMW 
WBITE  (6,360)  HESTLT 

SGJITO=SGUITO*XMS 

WBITE  (6  ,350)  SG JITC  ,  AO A  P 

FLUX=FLUX*XK2*2. 

TILLUM^TIILUM*  2. 

WBITE  (6,370)  F1UX,TILLUM 
IF  (NFLGSN.EQ.  1)  WHITE  (6,170) 

STOP 

===========-=  I/C  FORMAT  ST  AT  EM  EN 1 S  =======  =  =  ==== 

500  FORMAT  (50XcF20-0) 

510  FORMAT  (50A1,E20.0) 

520  FCBMAT  (50A1.I20) 

530  FORMAT  (6CA1) 

C 

130  FORMAT  ( 1  C  X, f  BLOCKING  ADAPTIVE  CPTICS  FACTOR 
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140 

150 

160 

170 

180 

185 

190 

195 

220 


230 


240 

250 

260 

270 

280 

290 

300 

310 


320 

330 


|  1 CX, 1 IRBADIANCE  AREA  AVERAGING  FACTOa  =  *,510.4,/) 

FORMAT  (1  OX,  '  TELESCOPE  DIMENSIONS: 1 ,/ , 

1  1CX, 'OUTER  DIAMETER  =  ',F10.2,'  CM',/, 

1CX, 'INNER  DIAMETER  =  ',F10.2,'  CM',/. 

10X, ' GAUSSIAN  WAIST  DIAMETER  =  ' , F 1 0 . 2 ,  *  CM',/) 

FORMAT  (1 CX, ' RMS  WAVES  DISTORTION  =  ’,£10.4,/, 

j  1CX,  'PRASE  CCRR.E1TATIO  N  LENGTH  =  ',£10.4  '  CM',/) 

FORMAT  (  1  OX ,  'ZENITH  LOG  AMPLITUDE  VARIANCE  =  ' , Z 1 2 . 4 ,/ ,  10X ,  ' REI A  '  , 
1  'LIVE  IRR  ADI ANCE  5EDU CT ICN  =  '  . E  12 . 4 , /} 

FORMAT  (  10X,  '  NCTE  =  =  =  — =  =  THE  ADAPTIVE  OPTICS  SYSTEM  USED  WILL  ', 

1  'PROVIDE  DEGRADED  PERFORMANCE  OVER  TEE  TILT  ONLY  CASE',/) 

1CI ,50  A1 , £20 .o) 

1CX,5CA1,I20) 

i1B1.9I.'1KP0T  DATA  FILE:',//) 

FORMAT  (1CX,90A1,/) 

FORMAT  (1H1.9X.  'PATH  ANALYSIS  RESULTS : ' , //, 

|  6  X  ,  S  T  r  F  ' ,  10X  .'RANGE'  . 9X,  ' SL  E  * '  ,6X, 'OMEGA' ,7X,  ' TIME'  , 

|  51,'ATHCS.  TRANSMISSION' , 10X,  '  R  H  G 1 , / 


FORMAT 

I  16 


1  IX,  ' FIOENCE  *  ;/, 


MAX  TREAD* 


4  X 


f  5X  ,  '  iXJE 


4 ’STEP'  5X  'DUE  10  '  ,  5  X  ,  '  AM  P  LOSS',  4X, 

'TEERMAI  *  ,4X  'SCALE'  5X,  MXW/CM2)  '  ,51, 

«  (KJ/CM2)  ',/.  1 0x  ,'  SPREADING'  ,  6 X.  'WITH  AO'  5X,  '  BLOOMING  '  ,//) 
FORMAT  (7X,F3.C,EX,PlC.3,5X,F3.4,5X,?6.2,5X,?6.2,oX, 

1  F7.5, 5X,F7.5 . 5X,F8.1 ,/) 

FORMAT (72,F3.C,5X,E10.:>,5X,Flt.5,5X,El0.5,5I,F10.4,5X,E10.4, 
i  5X.E1C.4,/) 

FORMAT  (  lOX, ' THIS  RON  -• ) 

FORMAT  (  12X,  ' INCLODES  THE  STATISTICAL  CICUD  MODEL') 

FORMAT  (  12X, ' INCLODES  THE  FULI  ZONAL  AC  MODEL  WITH  AN  ', 

1  '  ISOP1ANA 1ISK  MODEL') 

FORMAT  {  12X,  ' INCLUDES  THE  FULI  ZONAL  AC  MODEL  WITHOUT  AN  ’, 

1  'ISOPLAN ATISC  MODEL' ) 

FORMAT  {  12X,  '  INCLUDES  FULI  TILT  CORRECTION  KITH  TUT  ', 

1  'ISOPLANATISM') 

FORMAT  (121, 'IS  TEE  BASIC  CODE  KITH  ONLY  SOME  MEASURE  OF 

1  'TILT  CORRECT ION '  ,/, 1 2X, 'AND  WITHOUT  AN  I  SO PLA NAT IC '  , 

2  '  MODE I '  , /) 

FORMAT  (10X  'ZENITH  LOOK-AHEAD  ANGLE  =  *,E10.4,'  GRAD',/) 

FORMAT  (1E1,9X,  'CASE  DATA  AND  CALCULATED  FACTORS:',//, 


340 


350 


360 

370 


380 


1  GX, 
10X, 
1  GX, 
1CX, 
1  OX, 

1  cx, 
1  cx, 
1  cx, 

1  OX, 
10X, 
1  cx, 
1  cx, 

1  OX, 

FORMAT  (1 GX, 
1  2X, 
121, 
1  2X, 
1  21, 
122, 

F CBM AT  (1 CX, 
1  OX, 
10X, 

FORMAT  (1  OX, 
FORMAT  l/^, 

1  c  x# 

FORMAT  (lOlI 
1  2  X , 

i  1  2  X, 

,  1  2  X , 

I  12X, 

END 


TRANSMITTER  DIAMETER 
CAVITY  WAVELENGTH 
APERTURE  TOTAL  POWER 
TIMES  DIFFRACTION  IIMIT  = 

(WIDE  ANGLE  SCATTERING)  '  ,/, 
TRANSPUTER  ALTITUDE 
SATELLITE  ALTITUDE 
MAX  ZENITH  ANGLE 
FLIGHT  FATH  OFFSET 
RHCO 

WIND  VELOCITY 
HEIGHT  OF  GROUND 
OPT  SEEING  AT  550CA  =  . 

N  UMEEfi  CE  INTERVAL  STEPS  FOR:',/, 
ANGULAR  INTERVAL 
A  E  SCRFTICN  INTEGRATION 
RHC  CALCULATION 
MTF  CALCULATION 
THERMAL  BLOOMING 
2-SIGMA-F  BEAM  JITTER 
TURBULENCE  JITTER  E EJECT  ION 


, F10.2, '  CM' ,/, 
, F 1 0 . 6 , '  UM',/, 


, F10 .2, • 

,  F  1 0 . 2  ,  • 


M  ’  / 

KM  '  ' 


,  r  i  u  .  £.  ,  au  ,  /  , 

,  F  10 . 2 , '  DEG ' , /, 

, F  10 . 2 ,  '  KM ',/, 

, E 1 0 . 4 , '  CM't/. 

, F  10 . 2 , '  M/s£c(,/, 


,F  10 .2  , 

, E10.2, ' 


an C-S  lC ' 


r/) 


=  '  ,no,/, 

=  ' ,110,/, 

=  ',110,/, 

=  ',I10,/f 

=  '  ,  ?  1  0 . 2 , '  ORAD',/, 
=  1 ,210. 4,/, 


(RFEItOAL  =  INITAl  *  ADA  F) ' ,/) 

ISCPlAtiATIC  JITTER  (2-SIGMA-P)  =  ',£10.4, 
OX,' ========  PROPAGATION  SESOlTS  =  ===*='=* 


a  * - — - i 

INTEGRATED  FLUX  ON  TARGET 
T  Cl A  I  ILLUMINATION  TIME 
FOR  TEE  FULL  ZONAL  A-0  MCDEI: 
ZENITH  TRANS.  AT  *-0  SENSOR 
TARGET  RADIANT  INTENSITY 
ADAPTIVE  OPTICS  BANDWIDTH 
NONEEF  CF  ACTUATORS 


■/X 


=  '  ,  F 1 5. 6,  KJ/C  M2 ' ,/ , 
=  ' , F 1 5. 2,'  SEC') 

=/?,E10.4,/, 

=  ',£10.4,'  W/STER’,/, 

=  ',E10.4,'  HEP.T  Z  '  ,  /  , 

=  '  ,F10.0,/) 
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S  UEEOUTI N £  A  L  F  2  (A, 2) 

C 

C - CALCULATE  THE  TCTAI  MOLECULAR  ABSCEPTICN 

C 

DIMENSION  ALT  (2  0)  ,  AT  A  (20) 

C 

C******  THIS  TATA  FCR  CO  5? 9  FECK  AFWL  INFO,  REPORTEDLY  EASED  CN  AN 
(:**♦***  UPDATE  OF  TEE  MCCLATCHY  LINE  DATA.  THE  WAVELENGTH  FOR  THIS 
C******  TRANSITION  IS  4.99210  MICRCNS. 

C******  AT  MOSFHS  HE  IS  MI L-L AT ITU D t  SUMMER,  CLEAR  DAY. 

DATA  ALT/0. . 

i  n 

DATA  ATA/C. 

1  J  .  •+  V  Z.  _  U  ^  •  V  I  L.  —  U  J  ,  I.  /  U  J  ,  J  .  **  U  —  W  J  U  ~  U  - 

2  1.  061-05,6*  2  IE-06,  3. 61E-06,  C.0, 0*0,0. 0,0.0/ 

NI=  1 9 
I NDX= 1 
A  =  C . 


C 

C - 

c 


10 


IF  ((Z*3*279)  .GF.  1.E5)  RETURN 
H=Z*T. E-3 

SEARCH  FCR  THE  ALTITUDE  INDEX 


DO  10  1=1.  NL 
IF  (H.  LT.ALT (I) )  GC  TO 
I ND X=I 
CONTINUE 


10 


fi=  (H-ALT 

(INDX)  ) 

/  (alt  j: 

A  =  H*  {AT A 
A  =  A/  1 0  0  0 
RETURN 

(INDX  + 1 

)  -  AT  A  ( 

■AIT 


iJSSf" 


ESC 
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SUEROUTINE  ALFS  (Sf2) 

C 

C -  CALCULATE  THE  1CTAL  SCATTERING 

C 

DIMENSION  ALT  (20)  ,A1S  (20) 

C 

C******  TH I S  DATA  PC E  CO  5?9  FROM  A  FWL  INFO,  REPORTEDLY  EASED  ON  AH 
C******  OP  DA  I E  CF  TEE  KCCLATCHY  LINE  DATA.  T  BE  WAVELENGTH  FOR  THIS 
C  ***** #  TRANSITION  IS  4.592  10  KICRCNS. 

C******  ATMOSPHERE  1 5  MID-LATITUDE  SUMMER,  CLEAR  DAY- 

DATA  ALT/ 0.0 .  1-0, 1.0.3.  0.4. 0,5. 0.6.0,7.0.8.0,9.0,10.0,12.0,14.0, 
1  16.  6,  18.0.20.0,2 5.  0,30.0.35.0, 40.0/ 

DATA  ATS/  1. 4  1  E - 0 3 , 9  -  36E- 0 4,  1  .  C S E- J4 ,  5 . 3 7E- 05 , 4  . 5 4 E-  05 ,  3. 5 0 E- C 5 , 

1  2. 6 0E-05 ,1.792-0  5,  1  .  14E-05,  4.  5 9 E-05 , 3 . 6 4 E- 05 , 2. 9  IE- 05, 

2  2.62E-C5,2. coE- 0  5,2. 4  IE-05,  1 . 6 SE- 05 , 7 . 3 0 E-06 ,  1 . 06  E-05 , 

3  3-2  8E-06 , 1 .611-06/ 

N 1=  1 9 

C  -  Q  . 

IF  JfZ*3-279)  -GI.  1.E5)  RETURN 
E  =  Z*  T. 2-3 
INDX- 1 


C - SEARCH  F  CF.  THE  A  IT  IT  ODE 

C 

DC  10  1=1, NL 


INDEX 


10 


IE  (H. IT.iLT  (I) )  GC  TO  10 
INDX=I 


CONTINUE 
R=  (H-ALT 
5=a*  l ATS  , 
S=S/1000. 
RETURN 
INC 


(INDX)  )/  (ALT  (INDX +  1)  -  AIT  (IND  X)  ) 
(INDX+  1)  -  ATS  (INDX)  )  +A1S  (INDX) 
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nnn 


10 


SUBROUTINE  A3SCS3  (  N  ,  HT  ,  T ) 

CALCULI T £5  ZENITH  THE  TOTAL  INTEGRATED  MOLECULAR  TRANSMISSION 

COflflON  / ATflO/  EA I  EG 
DEI  H=( HA  T30- HT ) /N 
HEIGHT=Hl4DELH/2. 

A1FZ=0. 

DC  10  1=1, N 

CALL  ALPA  (ALF  ,  EE3G  HI) 

A LFZ= A  LP  Z  + AL  P 
BEIGHT=HEIGHT+DELE 
CONTINUE 
ALFZ=ALPZ*DELE 
T=EXP (-AIFZ) 

EETURN 

END 
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non 


SUBROUTINE  SCAT  (>,H1,T) 

- CALCULATES  ZENITH  THE  TOTAL  I NIEGRAT  EL  SCATTERED  TRANSMISSION 

COMMON  / ATMO/  HAT  EC 
DELH= (HA1M0-HT)/N 
HEIGHT I+D2LB/2 . 

AIPS*0. 

DO  10  1*1, N 

CALL  ALES  (ALS  ,  E  E I G  H 1 ) 

AIFS=ALPS+ALS 
HEIGHT=HEJGHT+E£LE 
10  CONTINUE 

A1PS=ALPS*DELE 

1=EX?{-A1ES) 

RETURN 

END 


non 


SUEROUTINI  RHC1EB  ( K , HT , H G , W  ,  £ E C) 

-  COMPUTES  YOEAS  2ENI1E  ATiiCSPEEEIC  COHERENCE  DIAMETER 

ccaaoN  /Also/  batcc 

TPI=4. ♦ARSIN ( U) 

CK=1PI/H 

CKSQ=CK*CK 

DEIH= (HAiaO-HT )yFICA1  (N) 

HEIGHT=HG*DELfi/^. 

SDM=0. 

DC  10  1=1, N 

CAIL  CN2E  (HEIGHT, CN2) 

SDfl=SUfl+CN2 
HEIGHT=HZ1GHT+CEIE 
10  CONTINUE 

CNI=SU  fl*  DELE 
E  E  N=2. 91 *CKSQ*CNL 
RHC=(6 .88/DEW)  **  (.6) 

EBC=RHO/2.1 

RETURN 

END 
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c 

c — 
c — 
c 


10 

c 

20 

C 

C — 

c 


40 

50 


SUBROUTINE  FARE  ID  (DX2,  F,  G,  M I,N4,HSATD, V A V Z, PTOTA L , PH A XO) 

COMPUTES  THE  F  AS  F J  El E  IRRADIANCE  OF  T  BE  APERTURE  DISTRIBUTION 
AND  THEN  COMPOTES  TEE  COR RESPC NDING  APEBTUEE  tttr 

DIMENSION  F  (HI)  ,G  (HI) 

COMMON  / AR8FL  D/  LI  A , LI  AC  E3, BEAC3Z,U0 
FI2=AB  SI N  (1.) 

F 1=2. *  PI  2 
T  WCPI  =  2 . *  FI 
C  K  =  T WO  PI/SAVE 
C  KR=CK/HS  ATD 
LIASQ=DI A**2 
DCESS0=D2A05S*  *2 

R FOCUS =1 . 2*HS A TL* WAV  E/S QE I (DI 2SQ-DOBSSQ) 

F  ACTOR  =  4 . 

F  ACTOR  =2 . 5 
F  ACTOR  =5 . 

RIIH=PACTCR*RFCCCS 
DX  0= (DIA-EIAOES) /2./ML 
DX1=RLIM/KL 
El=DXl/2. 

DC  20  1=1, ML 
SUM=0. 

R  0=DIA  OB  5/2.  +EX0/2. 

DC  10  J=1£ML 
Z=CKH*  fi  0  *  fi  1 
CALL  JO  J  2  ,  A  ) 

CALL  FIELD  (RC.OR) 

SOM=STJM+A*UB*BC 

B0=R0+DXC 

BY  DIVIDING  EY  FTCTAL,  THE  MTF  WILL  EE  UNITY  AT  THE  ORIGIN. 
F  (11  = ( SOM*CKB  *  EX  0)  ** 2/P  TOTAL 
fi 1  =  B 1+  DX  1 

COMPOTE  F  MAX  0  CN  AXIS 


SO  M=0 . 

B  0=DIAOBS/2. +DX0/2. 

DC  30  1=1, ML 
CALL  FIELD  (B  0  ,  OR ) 

S0M=SUM+  UB*R0 
B 0  =  E0+  DX  C 

P  B  AX0=  (SO  M*CK  F  ’♦DX  0)  *  *2 

C NORM-  FIXES  IFF  CTF  SO  THAT  TEE  MTF  INTEGRATION  PRODUCES  AN 
IEF1  VALUE. 

CNOBM=  (WAVE*HSATD) * * 2*P M A XO/P TCT AL 

CCBPOTI  TEF  KTf 

B2=D X2/2. 

DC  50  1=1, N4 
SDM=0. 

B  1  =  DXl/2. 

DC  40  J=  1 , ML 
Z=CKR*  E 1  *  £2 
CALL  JO  (Z ,  A) 

SOM=SO  M  +  A  *F  (J)  *E  1 
E  1  =B  1  +  DX  1 

G  m =SOM*1VOPI *DX  1/CNORfl 

R  5  =  B  2+  DX  2 

RETURN 

END 
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non  on 


SUBROUTINE  JO  <X,F) 

CCSPUTES  THE  ZERO  ORDER  BESSEI  FUNCTION. 


C 

c 


10 

20 


30 


40 


50 


60 

70 


DIKENS 10  N 
R  E  AL*3 
DATA 
1  6  36D' 


INSIOK  A  1  (16),A2(11),A3(11),A4(15),A5(11)  ,  A  6  ( 1  1 ) 
,*3  A  1  j  16)  ,  A$  (11)  ,A3  (1  1)  ,  A  4  1 15) ,  A  5  (1  1)  ,  Ao  f  11) 

A  1/-  1.0 D-  15.4.  l£- 14,-1 .  9  44D-  12,  7 . 6 4 8^ D- 1 1 , - 2 . c 

i-fl  -i  7Aicu^tr.qr-A  ^  9nAn*:3fnin-^  -n  t)0  62o1bo2 


679257D-9 ,7.6081 
bQ62o1bo206D-4, 4. 61  SI  800 


29869504/ 

DATA  A3/-1.0D-  15*4. CD- 15, -3.  3  T- 1 4 , 3.  0  ID-13,-3.  207 D-  12,  4. 220  ID -1 1 
17. 2719  2D- 10,  1.7S/245  7D-8,-7.4  1  44 9 &4 1 D-7 , 6 . §335  1 9942 6D- 5 , - . 03 1 1 1 1 ^0 
29210674/ 

DATA  A  4/1. ID-14, -5. 7  8D-  1  3 , 2  .  5  28  ID- 1  1  ,-9 . 4 242  1 D- 1  0 , 2 . 94 970  7  D- 8  ,  - 7 . 6 
1  1  75878  ID-7,1  - 58E7C1924D-5,-.  0 002 604 4 4  3 8  S3  4  9 , . 0 03 24 0 270  1 3 26 8 4. -. C 2 9 
2175 52 4806 1b4.  .  177709117239726, -.6614  43934134543, 1 .2879940986^768,- 
3  1.  1918  01  16054  122,  1.2967  1754  1  2  1C53/ 

DATA  A5/-1.0D-  15,5.  0D-1 5 , -4. 7L-  14,4.  7D- 13 ,-5.7  OSD- 12,6.8169D-11.-1 
1.87189  ID- 9, 6.  177  6  341-8,-3.9  87  28  430-6  ,.  0008 93 98 93330 8o , 2. 00  1  60 6  0c  17 
22003/ 

LATA  A6/1.0D- 15,-5. OD-15 , 3. 6 D -  14,-3.  26D-13, 3. 5 15D- 12/-4.G864D- 11.8 
1. 229  19  D-  1C,- 2.  C  5  597  6  ID-  8,  9.  1  3  6  6  1 5  2bD  -  7 , -9  - 1>277  23  5  4  9  26-  5  ,  .  0  9  3  5  55  5  ^4 
2139071/ 

ZERO  ORDER  EESSII  EUNCTION 
N  =  0 


IF  (N.  EQ.  1)  GC  TO  50 
IF  JABS (X) .GT. 6. 0)  GC  TO  20 

Y  = . 062  5*  X  *  X- 2 
E  =  0.  0 
B  F 1 =0 . 0 

DC  10  1=1,16  - 

EF2=BP  1 
EF  1  =  B 

B=Y*BP  1-EF2+A1  { I ) 

F=.5*  (B-EF2) 

RETURN 

Y=256/  (X*X)  -2 
AE=ABS  X) 

E  =  0. 0 
BP  1=0 . 0 
DC  30  1=1,11 
BF2=BP  1 
EF  1=  B 

B  =  Y*3P 1-EF2+ A2  (I) 

P=.5* (E-EF2) 

B  =  C.  0 
BF  1  =  0.  0 
DO  40  1=1,11 
EF  2=BP  1 
BF  1  =  B 

E  =  Y*BP  1-EF2+ A3  <I) 

Q  =  4*  (B-B  E  2)  /  A  E 

Y= AB- . 785398163357448 

F=.  797  884560 8C2865*  (F*COS  (Y)-£*SIN(Y))/SQRT(AB) 

IF  (ABS(X).GT.  8.0)  GC  TO  70 

Y=.0625*X*X-2 

B  =  0. 0 

BP  1  =  0. 0 

DC  60  1=1,15 

BF2=BP  1 

BF  1  =  B 

E=Y*BP  1-EF2+A4  (I) 

F  = . 0625*  (E-BP2) *X 
RETURN 

Y  =  256/  (X*X)-2 
A E=ABS  (X) 

B  =  0.  0 
EP  1  =  0 . 0 
DC  80  1=1,11 
E£2=BP  1 
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80 


90 


£  P  1  =  B 

B=Y*BP 1-EE2+ A5  (I) 
P=,5* (E-EP2) 

B  =  G.  0 
BP  1=0 . 0 
DC  90  1=1,11 
E  P  2=BP  1 


BP  1  =  B 

E  =  Y * B?  1- EP2 ♦ A  6  (I) 

Q  =  4*  (B-  3  P  2)  /A  E 
B=SIGN  (1 . C,X) 

Y= AB-2 . 35c 1 9  4  4  S  C  1 9235 


F=. 797884560 3 C2£ 65*  ( F*C  0  S  (Y)  -C*SIN(Y) J/SQBl  (AB) 
F  =  E*F 


EETURN 

END 
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noon 


SUBROUTINE  F I  x  1 


D  (B,CE) 

-  CCHPUTES  EXIT  APEFTUFE  FIELD  EISIHIB UII CH -  NOTE  IT  BUST  3E 

- A  XI-SY  E 1 BIC . 

CCHHOK  /ABBFL L/  EIA/CIAOBS/  BEAf'SZ/OO 
OB  =  0  • 

D=2. *H 

IE  (D.  GT .DIA.  CE.D.ll.CIAOES)  EETURN 
EflRAD=E£A£SZ/2. 

U  B  =  U0*  EX  £  (- (E/EBE&L) **2) 

BETUHN 

END 
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non 


SUEKOUIINZ  O0CS1  (E1,N) 

- COMPUTES  NORM  AIIZATI  Ch  CONSTANT  GIVING  ETOTAL  ENERGY  IN  BEAM 

COMMON  / ARBFLL /  IIA,riAOBS,3EAtSZ,UO 
PI  =  2*  *  AR  SIN  ( 1. ) 

EC=DIA/2 . 

BCSQ=RO**2 
RI=DIA  OB  5/2. 

BISQ=RI**2 
DX=ROSQ/N 
CST=PI *D  X 


u  u—  I  • 

S  Dft  =  0 . 

DC  10  1=1, N 
RCOTX=SQET  (X) 

II  {ROOTS. GT. EC.CE.ECOTX.  IT. HI)  GOTO  10 
CALI  FIELD  (RCCTX,UE) 

SUM=SUM+DE**2 
10  X=X+DX 

T  EMP=S  UM  *  CS7 
0  0  =  SQRT (ET/TEKE) 

BETURN 

END 


114 


non 


S  OE  ROD  TI  h  £  MTE ATM  (2  ,D,  H  HC,  A  D  AT  ,  f ) 

-  COMPOTES  ATMOSPEEHIC  MTF  FUNCTION 

11=1/ D 
A 1=1.-ADAF 
D  BHO=D/R  RC 

E  =  E  XP  (-XE**  (  1 . 66  67)  *  (1.  -  A  1*XD**  (.333  3)  )  *D£HO**  (1.6667)) 

RETURN 

END 
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non 


SUEBOOTIM  HI  I  JII  {  X  ,SI  G  J  IT  ,  V  t  \  E ,  F) 

-  COMPOTES  JITIEF  K TF  FUNCTION 

TPI=4.  * A  E  SIN  ( 1 .) 

CK=TPI/fiUE 

F  =  EXP (—  (CK*X*S1GJII)  *♦2/8.) 

feioen 

ENE 
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noon 


SDEHOUTIhE  MTFEC  <X,F) 

-  CCMPUTES  BAHUCfl  PHASE  MTF 

-  FUNCTION. 

COMMON  /IQ/  TEFLMI, VARBQ, 
TEMP=1 •-EXP(-  (i/SCALEQ)  ** 
F=  EXP (-V AF3Q*T IMF ) 

RETURN 

ENE 


.  ASSUMES  A  GAUSSIAN  CORRELATION 
KAYEEQ /SCALSQ 
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SDE ROUTINE  3Q I R  E 1  ( A  ,  N4  ,  D  X,  B  21) 
D I £ENS 10  K  A  (  N  4  ) 

SUK=0. 

X=DX/2 . 

DC  10  1=1, N4 
C AIL  MTFEC  (X,E) 

S0M=SUM+ A  (I)  *f  *X 
X  =  X  +  DX 

10  CONTINUE 

RE1=SUM*EX*6.26E1852 

RETURN 

END 
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SUEfiOUTINE  DE1UV  (A,DX,N4) 

C 

C - IF  THE  RHS  WAVES  CF  PHASE  DIST  CE7I0N 

C -  THEN  THIS  SOUTINE  WILL  DETER SINE  THE 

C -  PEODUCES  TEE  SPECIFIED  TDEIHT. 

C 

C  CflflON  /EQy  TEFlHIfVAEBQ, WAV  EEC/SCAL 
DIHENS  10  N  A(N4) 

R  EIO= 1 •/TDFLMI<*2 
V  AEEQ=-A ICG j  a  EX  0 ) 

DVEC=VAEEC/10. 

TEST=. OOS 
NSIGN0=  *1  . 

10  CONTINUE 

CALL  BQIEEL  (  A  ,  N  4  f  D  X  ,  HEL) 

DEE1=B El-EELO 

IF  (A3S  (LEELWEEIC.it. TEST)  GC  TO  40 
IF  (DREL-IT* J. )  GC  TC  20 
NSIGN 1  =  1 
GC  TO  30 
20  N  SIGN  1 =-  1 

30  CONTINUE 

IF  ( NS IG K  1. H E. NSIGN 0)  DV3Q=DVE£/2. 
VAEBQ=VAEBQ+DVEv*FSIGN1 
NSI GNO  =  NSIGN 1 
GC  TO  10 
40  CONTINUE 

T S0PI=  6 *  253 1  6  S 2 
WAVEBQ=SCET ( VAFBC)/l«OPI 
EETUBN 
END 


ABE  NCT  SP 
AF  EROPEI AT 

C 
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tot>J 


non 


SUBROUTINE  JIT  (A  ,DX  ,  N4  ,SIGJI1 ,  WAVE) 
COMPUTES  ARRAY  FCF  JITTER  MTF. 


DIMENSION  A { N  4 ) 

X=DX/2 . 

DC  10  1=  1  ,N4 
CALL  MTFJIT  (X,SIGJIT 
A  (I)  =F 
10  X  =  X  +  DX 
RETURN 
END 


SAVE,  F) 
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nnn 


SDEfiOOTINZ  BA  (A,IX,N4) 

- COMPOTES  ARRAY  FOB  E  EAM  QUALITY  MTF. 

DIMENSION  A  {  N  4  ) 

X=DX/2 - 
DC  10  I=1/N4 
CALL  MTFEQ  ( X  ^  P ) 

A  (I)=F 
10  X=X+DX 
EETOBN 
END 
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nno  onnnn 


SUEROUTINE  CN2H  (  £  E  I  GHT  ,  CN2 ) 

CACULATE  ATMCSEfcEEIC  VERTICAL  TURBULENCE 

- HUFNAGELE  LATEST  £CD  EL-GO  CD  0  MY  ABOVE  3  KILOMETERS 

CN 2  =  2.  2*  M.  E-5  ♦  10  .  +  *  (-.  3  )  *HEIGBT)  **  1  0.  *EXP  (-HEIGHT/  100  0-  ) 
CN2=CN  2  +  i  *2- 1 6  *  FXF  {-£  El  GHT/15CG.) 

CN2=CN2*2.7 

-  THIS  MODIFICATION  IS  USED  to  include  turbulence  at  lower  alts 

CN  2=CN  2  +  1.2-1  4  *ZXE  (-HEIGHT/!  00.) 

RETURN 

END 
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S  UBROU  71 N  E  AV  (N  , VO , EHIDCT, HT  , CC5WND ,0020x3, -) 

-  COMPUTES  TEAT  FART  OF  THERMAL  E1CCMING  TEAT  CHANGES  *ITH  SLANT  FATh 

COMMON  / A1MO/  EATfiC 
COSSQ=CC£SND 
DEI H= (HAT  SO- HI ) /N 
RG=DEL  H/ 2 • 

HEIGHT=HT*DELE/2. 

VOCOS=VO*COSSC 

SUM=0. 

SUM AS=0 . 

DC  10  1=1, N 

CALL  ALFA  {ALA, HEIGHT) 

CALL  A  LF  £  ?AL£ -HEIGHT) 

SUftAS=SUKA3+AL£+AIA 
ALCSS=EXF  (-SUM  A S* LEI  E/C  OS CHG 1 
SUM=SU  M  +  A LOSS  *  ALA/  ( V CCO S  +  HG * ? HI  DOT) 

RG=RG+DFIE 
HEIGHT=HFIGHT+E£LH 
10  CONTINUE 

SUM=SOM*LELE 

E=SUM**2 

RETURN 

END 
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c 

C' 

c 

c 


10 

20 


SUBROUTINE  PEI  FES  (EI,MLH ,W AVI  ,PH) 

COMPOTES  THAT  F  A  El  OF  THERMAL  BLOOMING  PHASE  THAT  DOES  NOT 
CN  SLANT  RANGE 


COMMON  / AR3FLC/  FI A , CIA  0 ES, BE ACSZ ,00 
DIMENSION  PH  (  M I  E ,  ML) 

HLM=ML-1 
MLHM=MLH-1 
DI ASQ^DI b**2 
BADSQ=DIASQ/4  . 

PI4=AR  SI N  ( l - ) / 2 . 

CK=8.*PI4/»AVE 

SC ALE=  CK  *  (2.72E-4*.4)/(l .01£5*1.4) 
RCSQ=(CI£/2. ) **2 
RISQ= (DIACBS/2.) **2 
DEIN=DIA/MLM 
X I=DEL  N/z . 

YI=-DIA/2. +DE1N/2  . 

X  =XI 

DC  20  I=1,MLHM 


Y- YI 
S  UM  =  0 • 
xsg=x*x 
DC  10  J=  1  ,  ML  M 
BSQ=Y* Y+XSQ 
R  =  SCRT  (RSC) 

E  ACTOR  =  0  . 

IF  {RSQ.  IT.RCSC.ANB.  rSQ.  GT.  RISC)  CALL  FIELD  (E 

F  ACTOR  =  F  ACTO  R  *  *  2 

SOM=SUH+FACTOE 

PH  II- J)  «SCfl*SC A1F*EEIN 

Y  =  Y  +  DElN 


CONTINUE 


X  =  X  +  DE  LN 


CONTINUE 

RETORN 

END 


FACTOR) 


DEPEND 
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AND  FOCUS 


C 

C— 
c — 
c 


c 

c 

c 

c 

c — 
c 


c 

c— 
c  — 
c 


10 

20 

c~ 


c 


c 

c 

c 

c 


30 

40 


SUBROUTINE  RESEHS  (El  ,MLH  ,PH ) 

CALCULATES  TEE  RESIDUAL.  PHASE  AFTER  TEE  MEAN, TILT, 
HAVE  BEEN  REMOVED. 

DIMENSION  PH{K1H,HL) 

COMMON  / ARBFLE/  Cl A ,DIAOES, BE AESZ ,U0 
MIHM=MLH-1 

DELN=D IA/ELM 
X I  =  DEL  N/ 2  . 

Y I=-DI A/ 2 . +D  El N/2  . 

ECSQ= (DIA/2. ) **2 
B  IS  Q= ( DI ACBS/2  .) **2 

COMPUTE  NOBEALIZAIION  CONSTANTS 

C  AIL  COE  FS  {AC,  Al,A2,A3,  A4,MLK,MLHM,  DELN) 

COHPUTE  EXPANSICN  COEFFICIENTS 


S  U  E 1-0 . 
S  UM2-0  • 
S  UM3  =  0  . 
SUM4=0. 


X  =  XI 

DO  20  1=  1 , HLH E 
X SQ=X*  *  2 
Y  =  YI 


DC  10  J  =  1 , HLH 
YSQ=Y**2 
BSC=ISC+  YSQ 

IF  (BSCj.G1.BCSC.OE.ESC.LT.EISC)  go  tc 

CALL  FIELD  ( S C B1  { E S C) , OR) 


10 


DETERMINE  EXPANSICN  COEFFICIENTS  RELATIVE  TO  A  UNIFORM 
WEIGHTING  FUNCTION. 


F  1=0. 

IF  JABS  (CB) . GT.  1.1-2) 
F  ACTOR  =  PE  (I, J)  *FT 
SUE  1= SUM  1  +  FACICB  *  X 
SUH2=SUfl2+FACTCR*1 
SUM3=SUM 3 +PACTCR* ESQ 
SUM4=SUM4+FACTCR 
Y  =  Y  +DE  LN 
X=X+DELN 


FT=  1. 


—  ADJUST  FOR  KALE  PLANE  INTEGRATION 
SUM  1=0 . 

S  UM2=S  UM  2 * 2. 

SUM3=SOH3*2. 

SUH4  =  S  UM  4*2 . 


DEINSQ=DE1N**2 
PHMEAN=AC*SUM4  *DE1NSC 
P  HTLTX=  A  1*SUH1*DELNSC 
PBTLTY=A2*SUM2*DE1NSC 
PBFOC=  (A3*SUM3*A4*SUM4)  *DELNSQ 


SUBTRACT  THE  MEAN  ,1111, AND  FOCUS  CURVATURES. 


X  =  XI 

DC  40  1=  1 ,  MLH  M 
XSC=X**2 
Y  =  YI 

DO  30  J=  1  ,MLfl 
YSQ= Y*  *  2 
BSC=XSC+ YSQ 

IE  (RSQ.G1.ROSC.OB.BSQ.lt. RISC)  GO  TC  30 

PH(I,J)=FB{I,J)-FETLTX*A1 *X -p BTLT Y* A  2* Y -F HFOC*  ( A  3  * RS Q* A4)  - P HH E A N * A 

Y= Y+D2LN 
X  =  X  +DELN 
RETURN 
INC 


125 


10 

20 
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SUBROUTINE  CCEES  ( AO , A1  ,  A  2,  A  3 , A4 , MLM , MIH K , DS LN) 

COMPUTES  TILT  AND  FOCUS  NORM  A 11ZTI0N  CONSTANTS 

COMMON  /ABBFLI/  El  A,  DIAO  £S,  B  E  AMSZ  ,  UO* 

EOSQ= (DI A/2. )  *  *2 
RISQ= { Cl  A CBS/ 2 . ) *  *2 
DEINSQ=BFIN**2 
YI--DI A/ Z.+DE1E/2 . 

X I=DEL  N/2 . 

X  =  XI 
SUM1=0. 

S  U  M  2=0  - 
S  UK3=0 . 

SUM4=0  - 

DO  20  1=1, ML H S 
XS£=X*  *2 
Y  =  YI 

DO  10  J—  1  ,  ML  M 
ISC=Y**2 
ES£=XS  C+  YSQ 
ESCSQ=  BS£**2 

IF  (RSQ. GT.RCSC.CB.ESQ-IT. RISC)  GO  TC  10 
CALL  FIELD  ( S C FT  (  IS C ) ,0 E ) 

-DETERMINE  NORMALIZATION  COEFFICIENTS  RELATIVE  TO 
-WEIGHTING  FUNCTION. 

FT=Q . 

IF  (ABS  (CE)  .GT  .  1.  E-2)  FT=1. 

FACTOR=FT 
SUM1=SUM 1+FACTCB 
SUM2=SUMZ*FACTCF*XSC 
SUM3=S  UM3  +  FACTCF*XSC 
SUM4=SUM<J  +FACICB*ES£SQ 
Y=Y  +DELN 
X  =  X  +  DE  LN 

-  ADJUST  FOR  EALF  PLANE  INTEGRATION 

SUM 1=S  UM  1*2. 

SUM2=S  UM  2*2. 

SUM3=SUM  3*2. 

SDM4=SUM4*2. 

A  0  =  SQR  T  (1./(SUC1*EELKSQ)  ) 

A  1  =  SQR  T  1./ (SUM2*DE1NSQ) 

A 2  =  SQB  T  1./ (SUK3*CEINSQ)  ) 

E  =- SUM  1/ (SUM2  +  SUK3) 

A3=SQB T  (E*E/ (  (E*B*5UM4-SUM1)  *Er INSQ)  ) 

A  4  =  AJ/B 

RETURN 

END 
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A  UNIFORM 


c 

c — 
c 

c- 

c- 

c 


10 

20 

30 

C- 


C 

C 

C 

C 


SUBROUTINE  PHVAS  IS  I  ,  MLH  ,  JTOT  AI ,  PH,  S  IGSQO) 
COMPUTES  THE  TfcERKAI  BLOOMING  EHASE  VARIANCE 


VARIANCE  IS  CCEPUTED  LIKE  THE  STEEHL  RATIO  ACCORDING  TO  THE 
FIELD  AS  A  WEIGHTING  FUNCTION. 


DIMENSION  P  H  (  P 


3Q=  DI 
SQ*jEI 
5  =  ML  - 1 


R  CSQ= 
EISQ 
MLfi  = 
MLHM=M1H 


/  m, 

lb,pl) 


A/2. )  **2 
AC3S/2 .)  **2 


DILN=DIA/PLH 

XI=DELN/2. 

Y I=-DI A/ 2  -  +  XI 
X  =  XI 

suai=o. 

SUE2=0 . 

SUM3=0. 

DC  30  I=1,MLHM 
Y  =  YI 


EIAO  E  S, BEAPSZ,U0 


XSQ=X*X 
DO  20  J=1,ML3 
BSC=XSQ+  Y*Y 

IF  (ESQ- GT. RCSC.CB. ESC- LT.BISC) 
CALL  FIELD  ( S  C  RT  (ESC) /UR) 

SUM  1=3  UM 1 +Ua 


TC 


SU32=SUM2+UR*FE 
SU33=SU32+U3*F t 
CONTINUE 
Y=Y+DELN 
CONTINUE 
X  =  X  +  DE IN 


I/J) 

I.J) 


**2 


CONTINUE 

—  ADJUST  FOB  KALE  PLANE  INTEGRATION 
SUM 1  =  SUM  1+2- *  D  El N  *  *  2 
SUK2=S  UM  2  +2.  *CEIN**2 
SUM3=S  DM 3  +  2. *  BEIN'**  2 


10 


NORMALIZE  HITH  RESPECT  TO  INTEGRA  I  OF  THE  FIELD. 

SIGSQ0  =  SUM2/SUK1-  (SUM3/S  UM 1 )  **2 

RETURN 

END 
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10 


S  UEROU  TINE  ISCTRB  (N ,  HT  , HG,  R ,  ISCANG) 

-  COMPUTE  2ENITK  ISCF1ANATIC  ANGIE  BASED  UPON  D . 

—  DEE  INI  TICK 

COMMON  / A I MO/  EATMC 
REAL  ISOANG 
TPI=4 . *  A  B  SIN  (1.) 

CK=TPI/H 

CKSQ=CK*CE 


CEIB=(HATMO-HI)/N 

HEIGHT=HG+DELH/2* 

D  BI=DELH/ 2 . 

S  UM  =  0 . 

CC  10  1=1, N 

CALL  CN2E  {HEIGHT eCN2) 

S0M=SUtt+CN2*DHT** 1.  6  1 66 6 6 67 

H El GHT= HEIGHT +EEIE 

C  H1  =  DHT  +  C  ELH 

CONTINUE 

C  N1=SU  M*  T  ELH 

DEN =2* 91 *CKS£*CNL 

ISO ANG= (6„88/EEN)**.e*.314 

RETURN 

ENC 


L.  FRIEDS 
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c 

c 

c 

c 


10 


20 

30 


40 


50 

60 


70 

80 


90 
100 
1  10 
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C 

C 

C 


130 

c 

c 


S  UBROU  T I h  £  I  S  C  £  I A  ( J C  ,X  N  U  ,R  ,  S  i  CC  M3,  W  A  VE  ,  X I  NT  ,  HT  ,  HG,TISC,  N4  ,11) 

THIS  SUBROUTINE  COMPUTES  THE  ISOFLANAIIC  MTF  FOE  A  FULL  AO 
■C CMP £N  SATIE  SYS  IE*  DEVELOPED  EE  DL  F  El  EE 


DI  KENS  ION  3(23)  ,C  (19)  ,C  2  (  19)  ,  C  4  (19) 
ISC  ( N 4 ) 


DIMENSIONS  j  19}^  X  MU  (20)  ,  C N 2  (23) 


DIMENSION  T 

COMMON  /  AT  rt<5/  ^  ii  A  T 
IF  ( J  J  .  G  E  -  1)  GC  TC  3  C 
N  B  =  zO 
N  =  1 8 

DELH= (HATflO-Kl )/NE 

HEIGHT =HG+DELE/2- 

DC  10  K= 1 , NH 

CALL  CN2K  (HEIGHT  , CN  SQ) 

CN2  (K)  =CNSQ 

H  (K)  =HEIGHT 

HEIGHT  =  HE IGHT+DELE 

PI=  3.  1415926 

A  13  =  1 . /3  . 

A 5 6=5.  /6 . 

A572=5./T2. 

A  1  43=  1  4.  /3- 
A9  127=91  ./27. 

A73=7./3. 

A 53=5. /3  . 

A  1 2=1.  /2 . 

D  PBI=P 1/2 ./N 
PHI  =  DPBI/2- 
DC  20  J=  1  ,  N 
)=CCS 


P^I=PHI+£FHI' 


J)=CiJl*C((J) 


(J) 


GO  TO 
GC  TC 


60 

80 


C 4  jj  =C2  (0) 

CONTINUE 
JJ=JJ+  1 

XA=XNU*SECOMG/F 
DO  40  J=  1  ,N 

d(!J)iio’k  =  i,hh 

X=H  (K)  *XA 
XX=CN2  (K) 

IF  (X. LT .0.1) 

IF  (X.GT.10.) 

X  2=X*X 
DO  50  J=1,N 

F=1.+X** A 53- A  1 2*  (1.+2.*X*C(J)  +X2)  **  A  56- A  12  * (1.-2.*X*C(J)+X2)**A56 

S  ( J)  =S  ( J  )  +  F*  X  X 

GO  TO  IOC 

X  53=X*  * A  53 

X2=X*X 

X  4  =  X2  *  X2 

DC  70  J=  1  .N 

F  =  X53- A56*X2* j  1.-A13*C2  ( J ) )  ♦ A572*X4*  ( 1 . -A  1  4 3 *C2  ( J )  ♦ A 9 1 27*C 4  ( J )  ) 

s  (j)  =s  m+F*xi 

GO  TO  100 
X 13  =  X* *  ( —  A  1 3 ) 

Xl 3=1 1 3/  (X*X) 

DC  90  J=1,N 

F=1.-A56*X13*  (  1.-A13*C2  (J))  +A572*X73*  (1  .-A143*C2  (J)  +  A9 127 *C 4  (J)  ) 

S  (J)  =S  (J)  +F*  X  X 
CONTINUE 
CONTINUE 


XMU  (1)  =1  ./XA 
DO  120  J=  1 ,  N 

XflO  (J  +  1)  =114.  8e*DELB*S(J) 
NOW  INTEGRATE  CVEE  PHI 


CSI=B**A53*SECC  MG/WAVE* *2 
SUM=0. 

DC  130  1=1, N 

SUM=SUM+EXP (-XM0  (1+1) *CST) 
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TISO  (II)  =  XINI 

FZIUEN 

END 
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nnon 


SUEBOUTINE  FNEEEC  (  Ej£  D,  E  HO ,  D I  f.  ,R  HOU  ,  D  X ,  N4 , 1 R  RM  TF) 

- *  THIS  SUBROUTINE  DETEBttINES  WHAT  VALUE  CF  EHO  WOULD  PRODUCE  THE  AC 

-  CCEBECIEL  STHEHI  AS  If  THERE  WIRE  NO  COMPENSATION  AT  ALL 

BEAL  IBBETF  (N4) 

B=BHOU 
D£- BHO  U/ 2 • 

N  SIGNO  =- 1 
NSIGN=-1 
10  CONTINUE 

CALL  TREREL  (IfRKTE, El A ,N4/DX/SE,R) 

If  {A8S  (SF-EXDJ/EXC.  IT.  .  1  )  GO  TO  40 
DIf  =  SR  — JtXI 

If  (DIF. II. 0. )  GO  TO  20 
N  S IG  N-=  1 
GO  TO  30 
20  NSIGN=-1 

30  CGNTINDE 

If  (DIF.GE.O.)  DE=DH/2 
N  SIGNO  =  N  SIGN 
B=B-DS*N£IGN 
GC  TO  10 
'  40  CONTINUE 

R  HO=B 
RETURN 
END 
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noon 


SUBROUTINE  THEE  El  { I RP.MT  F  ,DI  A  ,  N4  ,  DX  ,  IB  El  ,R  HO) 

- THIS  S  UBfiCUTI N  E  CCMFC1ES  THE  RELATIVE  INTENSITY  DUE  TO  A  TOTAI 

- TILT  OERBECTEI  SYSTEM  FOB  TUREULENCE  CNIY. 

BEAL  IERKT F(N4) 

P 1=  2.  *  AR SIN  {  1  .) 

S  0  M  =  0 . 

DX2=DX/2. 

X  =  DX2 

DC  10  1=  1 ,  N4 
F 1  =  IRR  MTF  (I) 

CALL  HTFATK  ( X , E I A , R HO , 0  -  ,  F  2 ) 

su«=s 

X  =  X  +  DX 

10  CONTINUE 

TREL=SUM*2.*FI*DX 

RETURN 

END 
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c 

c 

c 


c 

c 

c 

c 


S  U3B0U  TIM  TL1ISC  (D 

THIS  SUBROUTINE  CCHF 

COttHON  / ATMO/  EATHO 
N  =  100 
E56=5. /6 - 
Ec^=5  /3 

DS=  (HATMC-HTBA KS) /N 
H  =  HA3V  GD  +ES/2 
SUB=0. 


,OMX 

UTES 


OH  Y 
THE 


BESTLT,  SECCHG ,HTE  ANS,K ABVGD) 

BESIDU A  I  TILT  DUS  TO  ISOPLANATISK 


DC  10  1=1, N 
CALL  CN2E  (H,CN2) 

Z=SECOHG*H 

AHG1= (  (D* 2*0 MX) **2  +  (2*0 HY)* *2) **E56 
ABG2=  {  (D-Z*ONXj  **2  +  2*0  MY)  **2[**E56 
A  BG3= (Z**2* (OH  2 **2+CHY*  *2) )  **  E56 
AKG4=D  *  *  E  53 

AEG5= {  (Z*CHX)  **2  +  (D-Z*OHY)**2) **Z56 
ARG6=(  (Z*CHX)  **2+  (D  +  Z*OMY)**2)**E56 

SUH=SUH+  (4. *AEG4+4. * ARG 3- AR G 2 -AEG  5- A BG 1-ARG6) *CN2 

H=B+DS 

CONTINUE 

THIS  IS  TEE  RESULTANT  2-SIGMA-F  TILT  DUE  TO  ISOPLANATISH- 

R£STLT=SUrt*DS*SECCMG*2.  9 1/D**2*2. 

RETURN 

END 
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SUE ROUTINE  RZITCT  (T,TR,TM) 

-  BE  HAVE  FOUND  TEE  DETAILED  flip  CODE  TC  GIVE  RESULT 

- WHICH  LIE  BETWEEN  AN  ESS  rlODIi  rCR  ELCCRING  AND  A 

-  HULTIPLICATIVE  APPROACH.  THIS  SUBROUTINE  IS  AN  AT 

-  AT  BETTER  MATCHING  THE  RESULTS  CF  TEE  RTF  CODE  BY 

-  AVERAGING  THE  ESS  AND  THE  M UITIPLIC ATI VE  RESULTS. 


T=  (IH«*TH)/2. 

RETURN 

END 


EH  FT 
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SUBROUTINE  BLOCK  (D,£,S,T) 

C 

C - THESE  BLOCKING  IRE1  KCDELS  ARE  CURVE  FIIIED  TO  GUISMTF 

C - HUNS  FOE  SIGSQ'S  GREATER  THAN  ABOUT  1.2.  THE  APE2TUP.E 

C - DISTRIBUTIONS  USED  HAD  OBSCURAT1CNS  OF  .1  X  THE  OUTER 

C - El AKETER-  IF  EETTEE  IS  REQUIRED,  EEC  CftflEND  BASELINING 

C - TO  THE  GUTSHTF  CODE  AGAIN. 

SOETS=SQRT JS) 

TEG=- . 0870d+5QETS*2.S1485+. 172 3* S 
TRG=1./TBG 

TEU=1. 2877-SQBlS*2.€U91+4  .096C2*S 
TBU= 1 . /TEU 


RETURN 

END 
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SUBROUTINE  SCI  NT  ISIGXZ  ,  N  ,HTB  ANS  , HGR  ND,  i»AVE) 

-  THIS  ROUTINE  COMPUTES  THE  VARIANCE  OF  TEE  LOG  AMPLITUDE 

-  FOR  A  ZENITH  ANGLE.  THE  RESUI1S  WILT  ONLY  BE  USED 

-  IN  THE  EVENT  TEAT  FULL  AO  IS  UTILIZE!  IN  THE  RON. 

COMMON  / AT  MO/  EATMC 
CK=6. 28/SAVE 
HTOTAL=BATMO-HTEASS 
DB=HTOTAI/N 
H=HTRANS+EH/2. 

C  5  6=5- /6 . 

S0M=0. 

DC  10  I=1.N 
H  5  6=H*  *C  5  6 
CALL  C  N  2  E  (H,  C  N  2) 

SUM=SOM+CN2*B56 
H=  E  +  DH 

10  CONTINUE 

C K7 6=C  K*  *  (7. /6  -)  *EH*G-56 

SIGXZ=SUM*CK76 

BETORN 

END 
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SU EROUTINE  SIiilIS  (T AMP , S IGXZ  ,  5EC0MG ) 

- COMPUTE  TEE  LOG  AMPLITUDE  SCI  El ILLAT ICN  VARIANCE 

- FOR  OFF-ZENITH  CONDITIONS  AND  THE  RELATIVE  IRHADIAHC3 

-  REDUCTION  WHEN  FOIL  AO  IS  USED. 

TAflP=EXP  (-SIGXZ*SECCMG**  (11./ 6.)  ) 

RETUHN 

END 
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